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PADRÕES NATURAIS 


A utilização do microcomputador na reprodução 
de padrões da natureza. Famílias de curvas. En- 
voltórias. Focos e cúspides. Teoria das catástro- 
fes. Padrões formados por pontos. Simulação do 
fenômeno da ressonância. Caos. Dinâmica de po- 
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PLANO DA OBRA 


“INPUT” é uma obra editada em fascículos semanais, e cada 
conjunto de 15 fascículos compõe um volume. A capa para enca- 
dernação de cada volume estará à venda oportunamente. 


COMPLETE SUA COLEÇÃO 


Exemplares atrasados, até seis meses após o encerramento da coleção, pode- 
rão ser comprados, a preços atualizados, da seguinte forma: 1. PESSOAL 
MENTE — Por meio de seu jornaleiro ou dirigindo-se ao distribuidor local, 
cujo endereço poderá ser facilmente conseguido junto a qualquer jornaleiro 
de sua cidade. Em São Paulo, os endereços são: rua Brigadeiro Tobias, 773, 
Centro; avenida Industrial, 117, Santo André; e no Rio de Janeiro: avenida 
Mem de Sá, 191/193, Centro. 2. POR CARTA — Poderão ser solicitados exem- 
plares atrasados também por carta, que deve ser enviada para DINAP — Dis- 
tribuidora Nacional de Publicações — Números Atrasados — Estrada Velha 
de Osasco, 132, Jardim Teresa — CEP 06000 — Osasco — SP. Não envie pa- 
gamento antecipado. O atendimento será feito pelo reembolso postal e o pa- 
gamento, incluindo as despesas postais, deverá ser efetuado ao se retirar a en- 
comenda na agência do Correio. 3. POR TELEX — Utilize o nº (011) 33 670 
DNAP. 

Em Portugal, os pedidos devem ser feitos à Distribuidora Jardim de Publica- 
ções, Lda. — Qta. Pau Varais, Azinhaga de Fetais — 2 685, Camarate — Lis- 
boa; Apartado 57 — Telex 43 069 JARLIS P. 

Atenção: Após seis meses do encerramento da coleção, os pedidos serão aten- 
didos dependendo da disponibilidade do estoque. 

Obs.: Quando pedir livros, mencione sempre título e/ou autor da obra, além 
do número da edição. 
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Com rotinas gráficas muito simples, 
podemos criar uma variedade 
surpreendente de desenhos — todos 
eles estruturados a partir de padrões 
fornecidos pela natureza. 


A utilização do computador para de- 
senhar a trajetória de objetos em queda 
livre sob a ação da gravidade foi exami- 
nada nos artigos publicados nas páginas 
766 e 781. Com os programas apresen- 
tados, ilustramos alguns aspectos de um 
velho ramo da Física chamado dinâmi- 
ca e vimos como traçar parábolas, cir- 
culos e elipses. Neste artigo, aprendere- 
mos a desenhar outros tipos de curva 
usando técnicas semelhantes. Rotinas 
gráficas bastante simples se encarrega- 
rão dessa tarefa. 

Nos últimos anos, a dinâmica rece- 
beu um novo impulso enquanto campo 
de pesquisa. Um dos motivos para isso 
foi o reconhecimento de que seus prin- 
cípios matemáticos não se aplicam so- 
mente ao movimento dos corpos sob a 
ação de determinadas forças. Cientistas 
ópticos, interessados na refração da luz 
das estrelas na atmosfera, químicos in- 
dustriais que estudam o desencadeamen- 
to de reações e biólogos preocupados 
com o crescimento de populações ani- 
mais que competem entre si — todos 
têm recorrido aos conhecimentos acu- 
mulados por esse ramo da Física. Além 
disso, os computadores — permitindo 
a repetição rápida de operações elemen- 
tares — favoreceram a revelação de fun- 
ções de uma complexidade estrutural 
que os antigos métodos de cálculo se- 
quer sugeriam. Com fregiiência, um pa- 
drão ou curva tem que ser repetido inú- 
meras vezes antes que possamos notar 
a presença de uma estrutura subjacen- 
te. Os programas abaixo demonstram 
isso. 


é FAMÍLIAS DE CURVAS 


Algumas curvas, isoladamente sim- 
ples, criam padrões bastante complexos 
quando são reunidas, constituindo uma 
família. Veja um exemplo em que se uti- 
lizam parábolas. 
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5 LET A$="": FOR N=1l TO 64: 
LET AS=AS+" ": NEXT N 
10 BRIGHT 1: BORDER O: 
PAPER O: CLS 
20 LET NMAX=81 
30 LET DELT=.05 
40 LET SX=170/S0R 3: 
175 
50 FOR N=1 TO NMAX 
60 LET A=PI*(-1+2*N/NMAX) 
70 PLOT 128,80 
B0 FOR T=0 TO 3 STEP DELT 
90 LET X=T*COS A: LET Y=T*( 
SIN A-T/2) 
100 IF Y<=-.4 THEN GOTO 120 
110 LET DX=SX*X+128: LET DY=ST 
*7+80 
111 IF DX<O OR DX>255 OR DY<O 
OR DY>175 THEN GOTO 120 
115 DRAW DX-PEEK 23677,DY-PEEK 
23678 
117 PRINT AT 19,0:A5 
120 NEXT T 


INK 5: 


LET 517= 


o 
Re 
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130 NEXT N 


Lo 


10 COLOR 15,1, 1:5CREEN £Z 
NM=81:PI=4*ATN(1) 
DE=.05 
5X=160/50R(3) 
FOR N=1 TO NM 
A=PT*(-1+2*N/NM) 
DRAW"BM12/,1187 
FOR T=0 TO 3 STEP DE 
X=T*COS (A) :Y=TA*(SIN(A)-T/2) 

100 IF Y>-.4 THEN LINE-(SX*X+12 

7, 118-8Y*Y) ELSE T=3 

110 FEXT T 

120 NEXT WN 

130 GoTo 130 


E (3) 


D)  HGRZ2 
20 NM = 81:PI = 4 * 
30 DE = .05 
Sx = 160 / SOR 
50 FOR N = 1 TO NM 


18Y=230 


ATN (1) 


(3) :5SY = 180 


- = sê 


60 A=PI*(-1+2%*%N/ NM) 
70 HCOLOR= 3:X1 = 139:Y1 = 118 
BO FORT = O TO 3 STEP DE 

90X =Tx*x Cos (A)J:Y=TAIS 
IN (A) TA 2) 

100 IF Y > - .35 THEN HPLOI 
Xx1,%Y1] TO 98X * X + 139,118 - Si 
* 7:xl = SX * % + 139:71) = 118 
- SY * Y: GOTO 110 

105 T = 3 

110 NEXT T 

120 NEXT N 


LO PMODE 4,1:PCLS 1:SCREEN 1,0 
20 NM=BL:PI=4*ATN(L) 

30 DE=.05 

40 SX=160/S0QR(3): 
50 FOR N=1 TO NM 
60 A=PI* (-1+2*N/NM) 

70 DRAW"BM127,11B” 

BO FOR T=0 TO 3 STEP DE 

90 X=TAáCOS (A) :Y=T*(SIN(A)-T/Z) 
100 IF Y>-.4 THEN LINE-(SK*X+12 


57=230 
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7,118-SY*Y), PRESET ELSE T=3 
110 NEXT T 

120 NEXT N 

130 GOTO 130 


O programa simula as trajetórias de 
jatos de água lançados por um esguicho 
de jardim, ou, em uma interpretação 
mais moderna, de nêutrons lançados de 
um tubo fino ligado a um reator. Em 
ambos os casos, a família de curvas é 
constituída por todas as trajetórias pa- 
rabólicas que emergem no mesmo pon- 
to, em direções variadas mas com a mes- 
ma velocidade. O padrão formado por 
essa família é a curva exterior que to- 
das as curvas tocam. À curva exterior, 
chamada de envoltória, é também, nes- 
te caso, uma parábola. Em balística, 
usa-se o termo parábola limitante para 
uma curva desse tipo, pois ela define os 
limites da região que pode ser alcança- 
da por projéteis que tenham uma deter- 
minada velocidade inicial. 

É importante notar que a envoltória 
é uma propriedade da família de curvas 
como um todo, não tendo nenhum sig- 
nificado para uma trajetória em parti- 
cular. Trata-se de um bom exemplo de 
que o todo pode ser mais do que a sim- 
ples soma de suas partes. 

No programa, a equação das curvas 
aparece na linha 90. X é a distância ho- 
rizontal; Y, a distância vertical; T cor- 
responde ao tempo (começando por ze- 
ro no instante do lançamento) e cada 
trajetória tem uma direção determina- 
da por A, o ângulo ou inclinação do lan- 
çamento. O programa desenha NMAX 
ou NM inclinações; tente mudar o va- 
lor dessa variável na linha 20, 


FOCOS E CÚSPIDES 





Este programa mostra como também 
as linhas retas podem constituir famílias 
com envoltórias interessantes. 


10 BORDER O: 
CLS : LET N=2 
30 LET Y0=80/N/N 

40 PLOT 0,-6-Y0*2 

50 FOR X=0 TO 255 STEP 2 

60 LET Y=169-Y0-YO*COS (N*2* 
PI*X/255) 

70 DRAW X-PEEK 23677,175-Y- 
PEEK 23678 

BO LET XT=X+2*Y0*N*PI*Y/255* 
SIN (N*2*PI*X/255) 

85 LET Y1l=0: IF XT<O THEN 

LET XT=0: LET Y1l=Y+255*X/(2*N 
*PI*YO*SIN (N*2*PI*X/255)) 

B6 IF XT>255 THEN LET XT=255 
: LET Yl=Y-(255-X)*255/(2*N* 
PIXTOXSIN (N*2*XPI*X/255)) 


INK 7: PAPER O: 


90 DRAW XT-PEEK 23677,175-Yl- 
PEEK 23678 

100 PLOT X,175-Y 

110 NEXT X 

120 GOTO 120 


a 


10 COLOR 15,1,1:SCREEN 2 
20 N=2:PI=4*ATN(1) 

30 Y0=80/N/N 

40 DRAW"BMO, "+STRS (INT (186-Y0*2 
)) 

50 FOR X=0 TO 255 STEP 2 

60 Y=186-Y0-YO*COS (N*2*PI*X/255 
) 

70 LINE -(X,Y) 

80 XT=X+2XYOXNXPIXY/255*SIN (N*2 
*PI*X/255) 

85 Y1=0:IF XT<O THEN XT=0:Yl=Y+ 
255*X/ (24NXPIXYOASIN(NA2*PIXX/2 
55)) 

86 IF XT>255 THEN XT=255:Yl=Y-( 
255-X) *255/ (24N*PIXYOXSIN(N*2*P 
1*X/255)) 

90 LINE -(XT,Y1) 

100 DRAW"BM"+STRS (X) +”, "+STRS(I 
NT(Y)) 
110 NEXT 
120 GOTO 


Dao 


10  HGB2 

20 N = 2:PI = 4 * 

30 YO = BO / N / W 
40 HCOLOR= 3:XX = 0:YY = 186 - 
vO * 2 

S0 FOR X = O TO 279 STEP 2 

60 Y = 186 - YO - YO * ZCOS (N 

* 2*%*PI XX / 279) 
70 HPLOT XX,YY TO 
BO XT =X + 2 * YO 


120 


ATN (1) 


XY 
NEPAL RS 


/ 279. * SIN (N*2Z2*PIY*XHOo/ 
279) 
85 Yl = O: IF XT < O THEN XT = 
O:Y] =Y+2Z79*KX/ (GRANT 
I* Y0* SIN(N*2*PI*XK/. 
279)) 


86 IF XT > 279 THEN XT = 279:Y 


L=T-(Zo=R) PESA ta 
N* PI * YO * SINQ(NY*Z* BI 
* X/ 279)) 


90 HPLOT X,Y TO XT,Y1 
100 XX = X:YY = Y 
110 NEXT 


10 PMODE 4,1:PCLSI:SCREEN 1,0 
20 N=2:PI=4*ATN(1) 

30 Y0=B0/N/N 

40 DRAW"BMO,"+STRS (INT(186-T70*2 
)) 

50 FOR X=0 TO 255 STEP 2 

60 Y=186-Y0-YOXCOS (N*2*HPI*X/255 
) 

20 LINE -(X,Y),PSET 

BO XT=X+2*Y0*N*PI*Y/255*SIN(N*2 
RPIXX/255) 

85 Y1=0:IF XT<O THEN KT=D:YZ=Y+ 
255*X/ (2*NXPIXYOXSIN (N*2*PI*XX/2 
25)) 

B6 IF XT>255 THEN XT=255:Yl=Y-( 


| ! = | 
| 
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255-X) *255/(24N*PIXYO*SIN(N*2*P 
1*x/255)) 

90 LINE -(XT,Y0), PRESET 

100 DRAW"BM"+STRS(X)+”,"+STRS(I 
NT(Y)) 

110 NEXT 

120 GOTO 120 


O programa simula o reflexo de raios 
de luz sobre uma superfície ondulada — 
como o sol batendo nas águas de uma 
lagoa. Esta família é composta somen- 
te por linhas retas, que formam ângu- 
los retos com uma curva senóide. A en- 
voltória consiste na curva constituída 
pelos pontos focais. Estes, no caso dos 
raios, são mais intensos nos vales da se- 
nóide, onde a envoltória tem pontos bri- 
lhantes chamados cuspides. A existên- 
cia de cúspides pode ser prevista atra- 
vés de um novo ramo da Matemática de- 
nominado Teoria das Catástrofes (esse 
nome dramático vem da aplicação da 
teoria ao desabamento de pontes e em- 
borcação de navios). 

No programa, a senóide é especifica- 
da na linha 60 e os raios de luz, na 80. 
O número de vales da senóide é N (me- 
nos no Apple). Tente mudar o valor de 
N na linha 20 (recomendamos N=1). 

Senóides também podem constituir 
uma família de curvas. Veja o progra- 
ma apresentado a seguir. 


10 BORDER O: 
CLS 

20 LET QM=S5QR (2Z*LN (3)) 

30 FOR Q=-QM TO QM*1.001 STEP 
QM/30 

40 PLOT 0,75+0*75/0M 

50 FOR T=0 TO 3*PI STEP .2 

60 LET X=Q*COS (EXP (-0*0/2)* 
T) 

70 DRAW (T*240/3/PI)-PEEK 
23677,75+(X*75/0M) -PEEK 23678 
BO NEXT T 

90 NEXT Q 


Da 


10 COLOR 15,1,1:S5CREEN 2 

20 QM=SQR (2*L0G(3)) :PI=4*ATN(1) 
30 FOR 0=-QM TO QM*1.001 STEP Q 
M/30 

40 DRAW 
0/0QM)) 
50 FOR T=0 TO 3*PI STEP .2Z 

60 X=Q*COS (EXP (-0*0/2)*T) 

70 LINE-(T*240/3/P1,100-X*90/0QM 
) 

BO NEXT T 

90 NEXT Q 

100 GOTO 100 


Cato 


10  HGRZ 
20 QM = 
= q * 


PAPER O: INK 7: 


"BMO,"+STRS (INT (100-0*9 


SOR (2 * LOG (3)):PI 


ATN (4) 














- QM TO QM * 1.001 





30 FOR Q = 
STEP QM / 30 
40 HCOLOR= 3:X1 = O:Y1 = 
Q* 90 / QM 

50 FORT =0OTO 3.5 * PI STEP 
E. 

60X = Q * 
1 2) * 1) 
70 HPLOT X1,Y1 TO T * 240 / 3 


100 - 


cos ( EXP (€C-0*70 


/ PI,100 - X * 90 / QM:X1L = T* 
240 7 3 /4/ PI:YI = 100 - X * 90 
f QM 

BO NEXT T 


90 NEXT Q 


a ad 


[ - 
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10 PMODE 4,1:PCLS 1:SCREEN 1,0 

20 QM=50R(2*L06(3)) : PI=4*ATN(1) 
30 FOR Q=-QM TO QM*1.001 STEP Q 
M/30 

40 DRAW"BMO,"+STRS (INT(100-0*90 
*QM)) 

50 FOR T=0 TO 3J*PI STEP .Z 

60 X=Q*COS (EXP(-0*0/2)*T) 

70 LINE -(T*240/3/P1,100-K*90/0 
M) , PRESET 

BO NEXT T 

90 NEXT OQ 

100 GOTO 100 


Simulamos aqui raios de luz atraves- 
sando uma fibra óptica cujo indice de 
refração varia ao longo de sua largura. 
Em uma escala microscópica, poderia- 
mos estar representando também “'tra- 
jetórias"" de elétrons entre dois planos 
de átomos de um cristal colocado sob o 
canhão de um microscópio eletrônico. 
A família é composta por senóides que 
começam no canto esquerdo do video, 
paralelas umas às outras; o periodo de 
cada curva depende de sua posição ini- 
cial. Embora esta família seja bem di- 
ferente da anterior, a envoltória também 
tem focos e cúspides. 

No programa, as órbitas são especi- 


TILT 
RCC 
At 


| Raios de luz em uma fibra óptica. 











ficadas na linha 60, que contém a dis- 
tância X do eixo principal no tempo É. 
As senóides são definidas por Q. Como 
está, o programa desenha trinta senói- 
des; para obter um número maior ou 
menor, modifique a linha 30. 





PADRÕES DE PONTOS 


Juntar curvas em famílias não é a 
única maneira de obter padrões interes- 
santes. Um outro caminho é seguir uma 
órbita por um longo tempo, de modo 
que uma estrutura complexa possa se re- 
velar, ainda que a fórmula matemática 
da curva seja relativamente simples. Pa- 
ra desenhar esse tipo de curva, não de- 
vemos fazer um traçado contínuo, pois, 
após alguns segundos, teremos preenchi- 
do a tela com um emaranhado de cur- 
vas. O mais conveniente é traçar a órbi- 
ta por meio de pontos, dispostos na te- 
la a intervalos regulares — como se a 
trajetória fosse observada sob luz estro- 
boscópica. Os três programas seguintes 
exemplificam essa técnica. A posição 
dos pontos na tela nem sempre equiva- 
le à posição real do objeto de estudo; o 
que o programa fornece, muitas vezes, 
é uma representação abstrata, na qual 
a coordenada horizontal corresponde à 
posição real e a coordenada vertical, à 
velocidade. 

Este programa leva vários minutos 
para completar o desenho: 


10 BORDER O: PAPER O: INK 7: 
CLS 

30 LET A=76.11 

40 LET ALF=A*PI/180: 
cos (ALF) 

50 LET S=SIN (ALF) 


LET C= 





| || | | 


NMAX=200 

M=52 

J=1 TO M 

X=0: LET Y=J/M 
N=1 TO NMAX 


60 LET 
O LET 
80 FOR 
90 LET 
100 FOR 
110 LET W=k 
120 LET X=X*C-(T-MX2M)FS: 
=Wí5+(Y-Wí*W) *C 
130 IF ABS (X)>4 OR ABS (7)24 
THEN GOTO B60 
135 IF X>1 OR Y>21 THEN 
150 
140 PLOT X*128+12B,Y*B5+85 
150 NEXT N 
160 NEXT qd 


LET Y 


GOTO 


ari 
q | , 


10 COLOR 15,1,1:SCREEN 2 
20 A=76.11:PI=4*ATN(1) 
30 AL=A*PI/180:C=COS(AL) 
40 S=SIN(AL) 


50 NM=200 

60 M=52 

70 FOR J=1 TOM 

BO X=0:Y=J/M 

90 FOR N=1 TO NM 

100 W=X 

110 K=X*C-(Y-X*X)*5:Y=WAS+(T-W* 
W) *C 


120 IF ABS(X)>4 OR ABS(Y)2>24 THE 
N 160 

125 IF ABS(Y)>1 OR ABS(X)2>21 THE 
N 140 

130 PSET(128+X*128,96-T7*96) 

140 NEXT N 

150 NEXT q 

160 GoTo 160 


———- — 
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10  HGRZ 


20 A = 76.11:PI = 4 * ATN (1) 
30 AL = A * PI / 180:€C = COS 1 
AL) 

dO S = SIN (AL) 

50 NM = 200 

60 M = 52 

70 FOR J = 1 TOM 





“Tlhas de pontos” ou linhas de força. 








a ii | " e ma =“ , pio ma = e” * 


BOX =0:7 = J/M 

90 FOR N = 1 TO NM 

100 W = X 

MO K=KSC=-(]7=-XAA "8 
TUAS +(T- WA) *C 
120 IF ABS (X) > 4 OR ABS (7 
) > 4 THEN 160 

125 IF ABS (X) > 1 OR ABS (7 


) > 1 THEN 140 


130 HPLOT 128 + X * 128,90 - 1 
* 96 

140 NEXT N 

150 NEXT q 


160 GOTO 160 


] 


10 PMODE 4,1:PCLSI:SCREEN 1,0 
20 A=76.11L:PI=4*ATN(1) 

30 AL=A*PI/180:C=COS (AL) 

40 S=5IN(AL) 

50 NM=200 

60 M=52 

70 FOR J=1l TO M 

BO X=D:Y=J/M 

90 FOR N=1 TO NM 

100 W=X 

110 X=K*C-(T-K*X) *5:Y=WA*5+(T-WA* 
w) FC 

120 IF ABS(X)>4 OR ABS(Y)24 THE 
N 160 

125 IF ABS(Y)>1 OR ABS(AX)>21 THE 
N 140 

130 PRESET(128+X*128,96-Y*96) 
140 NEXT N 

150 NEXT J 

160 GOTO 160 


O programa simula o movimento de 
partículas subatômicas (como prótons e 
elétrons) em um acelerador ou, ainda, 
linhas de força de um campo magnéti- 
co. Sua atuação, em termos matemaáti- 
cos, consiste em tomar uma série de 
pontos iniciais e movê-los pela tela, se- 
guindo esta regra: “faça uma rotação 
em torno do centro da tela, segundo um 
ângulo fixo, só que com uma ligeira mo- 





Uma regra simples produz o caos. 








dificação”. Sem esta “ligeira modifica- 
ção", todas as curvas seriam círculos 
concêntricos — e, de fato, as mais pró- 
ximas do centro se assemelham a circu- 
los. O efeito dessa modificação tem 
maior impacto sobre as curvas externas, 
que acabam formando um conjunto de 
“ilhas”. Se pudéssemos ampliar mais o 
desenho, distinguiriamos várias ilhas 
menores, distribuidas entre as que po- 
demos observar agora. 

No programa, o número de pontos 
iniciais é M, na linha 60: caso você não 
queira esperar tanto para ver o desenho 
completo, diminua esse valor. Cada 
ponto é desenhado NMAX ou NM ve- 
zes; O valor é especificado na linha 50. 
A regra fica nas linhas 100 e 110, que 
determinam de que modo as coordena- 
das vertical e horizontal serão modifi- 
cadas a cada repetição. O ângulo fixo 
de rotação é dado por A em graus, na 
linha 20; experimente valores diferentes 
(recomendamos 90). 

O programa apresentado a seguir ge- 
ra um padrão surpreendente, a partir de 
um ponto inicial: 


20 LET K=51.3: BORDER O: INK 
7: PAPER O: CLS 

30 LET X=1/PI: LET P=0 

40 LET A=1/BQR 5 

50 FOR N=1 TO 10000 

60 LET Y=X-.5: LET X=X+A-INT 
(X+A) 

70 LET P=P-Y 

BO PLOT X*255, P*K+B5 

90 NEXT N 


a f j 
k ! 
| 
F is 
LI. 


10 COLOR 15,1,1:SCREEN 2 
20 K=60:X=1/(4*ATN(1)):P=0 
30 A=1/SQR(S) 


) 





40 FOR N=1 TO 10000 

50 Y=X-.5:X=X+A-INT(IX+A) 
60 P=P-Y 

70 PSET(X*255,128-P*K) 
BO NEXT N 

90 GOTO 90 


Cao 


10 HGB2 :K = 40 

20 X = 1 / (4 * ATN (1)):P = 0 
30 À = 1) / SQR (5) 

40 FOR N = 1 TO 10000 

BO Ya N= .B:X = + A - 
(X + A) 

60 P =P -Y 
70 HCOLOR= 3: 
O - P*K 

BO NEXT N 


INT 


HPLOT X * 279,10 


10 PMODE 4,1:PCLS1:SCREEN 1,0:K 
=50 

20 K=1/(4*ATN(1)):P=0 

30 A=1/SQRIS) 

40 FOR N=1 TO 10000 

50 Y=X-.5b:X=X+A-INT(K+HA) 

OU P=P-Y 

70 PRESET (X*255,128-P*K) 

80 NEXT N 

90 GOTO 90 


Esse programa é uma simulação da 
ressonância, que ocorre quando dois fe- 
nômenos físicos têm frequências coin- 
cidentes ou múltiplas. Por exemplo, um 
dos fenômenos pode ser a órbita de um 
asteróide em torno do sol; o segundo, 
a perturbação dessa órbita pela atração 
gravitacional de Júpiter. Podemos nos 
perguntar: as perturbações causadas pe- 
lo grande planeta serão capazes de ar- 
rancar o asteróide de sua órbita? Tudo 
depende da razão entre as duas frequên- 
cias (a razão é obtida dividindo-se uma 
pela outra). As partículas dos anéis de 





Mudanças em uma população de peixes. 
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Saturno estão sujeitas a esse tipo de 
efeito. 

No programa, a razão é chamada de 
A e seu valor (menor que 1) está na h- 
nha 30. O número de repetições da res- 
sonância é 10000 (final da linha 40); re- 
duza esse valor caso não queira esperar 
muito. A transformação da ressonância, 
nas linhas 50 e 60, é obtida por meio de 
uma regra que modifica as posições ho- 
rizontal (X) e vertical (P) dos pontos na 
tela. 

Um delicado padrão em “'cortina”” é 
obtido com A = 1/SQR(5), ou seja, 
(0.4472136. Como esse valor não corres- 
ponde à razão entre dois números intei- 
ros (em termos do computador, trata- 
se de um número irracional), ele é não- 
ressonante. Para chegar à ressonância, 
temos que atribuir a A um valor que se- 
ja a razão entre dois números inteiros 
— como 9/20, que vale 0.45. Tente tam- 
bém um outro número irracional, como 
1/PI. Se quiser reduzir a escala vertical, 
diminua o valor de K. 


| CAOS | 





Observe que certas regras simples não 
produzem nenhum padrão. O próximo 
programa ilustra o fenômeno. +, 


10 BORDER O: PAPER O: INK 7: 
20 LET X=1/PI 

30 LET Y=1/PI 

FOR M=1 TO 10000 

SO LET X=X+Y-INT (X+Y) 

LET Y=K+Y-INT (X+T) 

70 PLOT X*255,Y*4175 

NEXT M 


a 


COLOR 15,1,1:SCREEN 2 
20 PI=4*ATN(1) :X=1/PI 
Y=1/PI 

FOR M=1 TO 10000 


50 X=X+Y-INT(X+Y) 

60 Y=X+Y-INT(X+T) 

70 PSET(X*255,192-Y*191) 

80 NEXT M 

90 GOTO 90 

Le) Lo) 

10  HGR2 

20 PI = 4 * ATN (1):X = 1 / PI 
30 7 = 1/ PI 

40 FORM = 1 TO 10000 

50 X=X+ Y- INTA(X+õgT) 

60 T=X+Y- INTÃ(XA+õY) 

70 HCOLOR= 3: HPLOT X * 279,19 


2: Wo M RDS 
80 NEXT M 


10 PMODE 4,1: PCLS1l:SCREEN 1,0 
20 PI=4HATN(1):X=1/PI 

sy Pe /PI 

FOR M=1 TO 10000 

50 K=X+Y-INT(X+T) 
Y=X+Y-INT(X+Y) 

70 PRESET(X*255,192-Y*191) 
NEXT M 

GOTO 90 


Esse programa simula o vaivém errá- 
tico das bolas de um fliperama, o mo- 
vimento das moléculas de um gás, ou 
qualquer outro sistema de órbita tão im- 
previsível que se torna impossível distin- 
gui-la de um processo aleatório. Con- 
vém notar, porém, que os 10.000 pon- 
tos desenhados na tela não são aleató- 
rios, resultando da repetição de uma re- 
gra determinística, ou seja, que não 
abriga nenhum elemento de acaso. A re- 
gra considera a tela como um quadra- 

o unitário onde as coordenadas verti- 
cais e horizontais variam entre zero e 
um. Ela opera em três etapas, nas linhas 
50 e 60 do programa. Primeiro, as coor- 
denadas X e Y do último ponto são so- 
madas para produzir um novo X; de- 
pois, esse X é somado a Y, resultando 
em outro Y; finalmente, se X ou Y fo- 
rem maiores que 1, um número inteiro 
apropriado é subtraído para que as 
coordenadas permaneçam no intervalo 
válido. 





POPULAÇÕES DE PEIXES 


Como a população de peixes de um 
lago varia ao longo das gerações? Isto 
depende dos fatores que condicionam as 
mudanças populacionais de uma gera- 
ção para a outra. Para definir uma re- 
gra, devemos considerar tanto a tendên- 
cia ao crescimento populacional resul- 
tante da reprodução, como a tendência 
à redução, determinada pela quantida- 
de de alimento disponível no lago. Con- 
forme o equilíbrio estabelecido entre es- 
ses dois fatores, a população pode ficar 
estável, alternar entre um ou mais valo- 
res, ou variar de tamanho ao acaso, ao 
longo das gerações. 

O programa dado a seguir ilustra O 
fenômeno. A posição horizontal, A, 
corresponde à razão entre disponibilida- 
de de alimento e taxa de natalidade, re- 
presentando assim as sucessivas gera- 
ções. À posição vertical Y corresponde 
ao tamanho da população, que aumen- 
ta ou diminui alguns pontos conforme 
as gerações se sucedem. O tamanho da 
população Y só é colocado na tela quan- 
do se encontra em equilíbrio, isto é, 
quando assume um valor estável ou al- 
terna entre dois ou mais valores, que são 





denominados pontos de atração. No 
Spectrum e no TRS-Color, um som de 
frequência proporcional ao tamanho da 
população é emitido. 


10 BORDER O: 
CLS 

20 LET 5=.03: 
LET NMAX=80 
FOR A=2.8 TO 4 STEP 5 

LET Y=1/PI 

FOR N=1 TO NMAX 

LET Y=A*YA(1-Y) 

IF N>NMIN THEN PLOT 255*( 
A-2.8)/1.2,Y*175 

BO SOUND .0075,Y*20 

90 NEXT N 

100 NEXT A 


io 


COLOR 15,1,1:SCREEN 2 
8=1/127:NN=50:NX=80 

FOR A=2.8 TO 4 STEP 5 
Y=.25/ATN(1) 

FOR N=1 TO NX 

Y=AFYAX(I-T) 

70 IF N>NN THEN PSET(255*(A-2.8 
)/1.2,192-1*191) 


INK 7: PAPER O: 


LET NMIN=50: 


"90 NEXT N 


100 NEXT A 
110 GOTO 110 


Dao 

10  HGRZ 

20 8 = 1 / 127:NN = 50:NX = 80 

30 FOR A = 2.8 TO 4 STEP 5 

40 7 = ..25 / ATN (1) 

50 FORN = 1 TO NkK 

60 YT=ART*H4I- 7) 

70 IF N> NN THEN HPLOT 279 * 
à — 2.8) É 1.2,6192 + FT. 7% 13% 

90 NEXT N 

100 NEXT A 


PMODE 4,1:PCLS1I:SCREEN 1,0 
8=1/127:NN=50:NX=B0 

FOR A=2.8 TO 4 STEP 5 
Y=,.25/ATN(1) 

FOR N=1 TO NX 

Y=AXYH(1-T) 

70 IF N>NN THEN PRESET (255* (A- 
2.8)/1.2,192-Y*191) 

B0 SOUND 1+255*Y,1 

90 NEXT N 

100 NEXT A 

110 GOTO 110 


'* Como você pode observar, os pontos 
de atração modificam-se drasticamente 
à medida que a disponibilidade de ali- 
mento aumenta. No início (lado esquer- 
do da tela), só há um valor de atração, 
indicando uma população estável. Isso 
significa que o alimento é suficiente pa- 
ra os peixes sobreviverem e se reprodu- 
zirem. Se a população crescer demais, 





haverá menos alimento e alguns peixes 
morrerão. Com a maior disponibilida- 
de de alimento para os sobreviventes, a 
população volta a crescer, podendo, em 
seguida, passar por uma nova fase de es- 
tabilidade. 

Subitamente, quando À assume um 
certo valor, a quantidade de alimento ul- 
trapassa determinado limite. Passam a 
existir, então, dois pontos de atração e 
a curva se bifurca — indicando que o 
tamanho da população agora se alterna 
entre dois valores. Mais tarde, quando 
a quantidade de alimento se torna ain- 
da maior, a curva se bifurça mais uma 


vez, aumentando o número de pontos 
de atração. Mais e mais divisões vão 
ocorrendo em uma sequência infinita, 
da qual o programa só tem resolução 
para mostrar os primeiros passos. As bi- 
furcações se acumulam e infinitos pon- 
tos de atração passam a corresponder a 
um certo valor finito de A. Em conse- 
quência, a população varia sem nunca 
se estabilizar. 

A descoberta dessa “árvore de bifur- 
cações”, que acaba num completo caos, 
despertou grande interesse na comuni- 
dade científica, por ter diversas aplica- 
ções. Ela retrata, por exemplo, a modi- 


ficação do fluxo de um líquido, que pas- 
sa de um estágio moderado (ponto de 
atração estável) para um turbulento 
(caos), conforme a velocidade aumen- 
ta. O mesmo efeito pode ser observado 
na fumaça de um cigarro, que sobe sua- 
vemente e de repente desenvolve espi- 
rais, terminando em turbulência. 

No programa, a lei de evolução está 
na linha 60. Para ampliar a resolução, 
altere a linha 20 de forma que S = .005, 
NMIM ou NM = 200€e NMAX ou NX 
= 300. No Spectrum e no TRS-Color, 
apague a linha 80 para aumentar a ve- 
locidade do programa. 
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AVAL/ span Nosso personagem não pode mesmo ; 
dae confiar em sua sorte. Surpreendido por 


Na uma avalanche tão logo começou a 
OS SALTOS DE WILUE iam 
= | = para escapar das pedras: saltar sobre elas. 


A parte do programa aqui apresen- avalanche. Vamos fazê-lo saltar, de mo- 410 REM ld de,259 
tada permite que Willie se movimente do que as pedras passem rolando sob 420 REM call 58970 
para cima e para baixo, no mesmo lu- suas pernas, sem atingi-lo. 430 REM ld de, 32 
gar — ou seja, nosso pobre personagem 440 REM add hl,de 


aprende, a tempo, a saltar sobre as pe- 450 REM 1d (57332) ,hi 
dras que-se aproximam. “E 




















460 REM ret 
4/0 REM m3d cp 4 
















































Willie já pode andar, mas isso não o 480 REM jr nz.mfj 
ajudará muito na hora de enfrentar a REM org 59336 490 REM Id a,0 
20 REM jmp ld de,6 500 REM ld (57335),a | 
30 REM ld h1,759 510 REM ld h1, (57332) 





40 REM call 949 520 REM ld de,32 

50 REM ld a, (57335) 430 REM abc hl,de 

60 REM cpl 540 REM ld bc,15616 

70 REM jr nz,mjb 550 REM ld a,45 

80 REM inc a . 560 REM call 58217 

90 REM ld (57335),a 570 REM j5yp 59153 

100 BEM 1d h1, (57332) SBO REM mfl ret 

110 REM 1d de,32 

120 REM sbc hl,de Para montar essa rotina na mes- 


130 REM ld (57332),hl mm 
140 REM 1d bc,57048 


150 REM ld a,40 


ma sequência da que demos no ar- 
tigo anterior, elimine o ret que foi 


160 REM Id de,259 colocado temporariamente como 
170 REM call 58970 último comando. Em seu lugar, en- 
180 REM ret tra o rótulo jmp, que significa salto. 
oº 190 REM mjb cp 2 Assim que começa a ser executa- 


200 REM 3r nz,mjc da, nossa rotina faz o computador 
REM inc a emitir um sinal sonoro. Para isso, 
REM ld (57335) ,a chama a rotina BEEP da ROM, que 
REM 19 bl, (3/3352) fica no endereço 949. Os parâme- 
tros adequados são mostrados nos 
pares de registros HL e DE. 


| REM ld bc,57000 
REM ld a,40 

REM ld de, 258 
REM call 58970 


REM ld de, 64 - 
REM add hl,de SUBINDO 





REM ld a,45 

REM ld bc,15616 A rotina dassalto compõe-se de 
REM call 58217 quatro partes; Gada uma delas é exe- 
REM ret cutada a seu tempo, quando a roti- 
REM mjc cp 3 na de movimentação é chamada. 
Sa RR ERG No último Artigo, mostramos 


REM Id (57335) ,a que o processador é direcionado pa- 
380 REM ld h1, (57332) ra a rotina Jump sé o conteúdo da po- 
* 390 REMÍId bc,57048 sição de memória 57335 não é ze- 
400 REM ld a,40 ro. Expliçãmos também como exa- 
minar qieciado para detectar qual- 
quer pressão sobreas teclas Me N. 





JAVA 


———= 
; m AS QUATRO PARTES E AJUSTAMENTO 
os DA ROTINA DE SALTO DOS APONTADORES DE TELA 
E * SALTOVERTICAL E WILLIE ESTÁ NO AR 
E SALTO DIAGONAL E DE VOLTA AO SOLO 
Es EFEITOS SONOROS E FIM DO SALTO 


) 


Se apenas a tecla N foi pressionada 
e Willie deve saltar, a posição de memó- 
ria 57335 é carregada com 1, que irá di- 
recionar o processador para a rotina 
aqui apresentada na próxima vez que o 
programa principal — que controla o jo- 
go — chamar a rotina de movimenta- 
ção. Caso as teclas M e N tenham sido 
pres onnaas, Willie deve saltar para a 

rente (em diagonal); 129 é, então, car- 
regado em 57335. Esta situação será de- 
vidamente analisada no próximo artigo 
da série Avalanche. 

Por hora, imagine apenas que a roti- 
na de movimentação de Willie consta- 
tou que houve uma pressão sobre a te- 
cla N, na última vez que foi chamada. 
Agora ela voltou a ser chamada, o que 
fez o processador direcionar-se exata- 
mente para nossa rotina. Depois da 
emissão do som que caracteriza o início 
de sua execução... 

O conteúdo da posição de memória 
57335 é carregado no acumulador ou- 
tra vez. Embora ele já estivesse ali quan- 
do esta rotina foi chamada, é provável 
que a rotina BEEP tenha utilizado o 
acumulador para realizar o efeito sono- 
ro do salto. Assim, convém carregar o 
valor em A novamente. 

Depois, o conteúdo de A é compara- 
do com 1. Naturalmente, esse valor es- 
tá presente; em consegiência, a próxi- 
ma instrução, jr nz,mjb, não é executa- 
da e o processador continua com inc a. 
Isto incrementa o conteúdo do acumu- 
lador para 2, valor que é colocado de 
volta em 57335 por Id (57335),a. Incre- 
menta-se, assim, o contador de salto, 
para que a rotina de movimentação 
acione a rotina jmp na próxima vez que 
for chamada. Nesse ponto, porém, a 
primeira parte do salto será deixada de 
lado; o processador executa a segunda 
parte. Como você deve ter reparado, ca- 
be à instrução cp 1, seguida da instru- 
ção jr nz,mjb, fazer o processador sal- 
tar para a segunda parte. 

A posição de Willie na tela é en- 
tão carregada dos endereços 57332 
e 57333 para o par de registros HL. 
O valor 32 é subtraído dessa posi- 
ção e colocado em DE. Com a exe- 
cução de sbc hl, de, o apontador se 
move para a posição imediatamen- 
» acima de Willie na tela. O resul- 
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tado é carregado de volta nos endereços 
$7332 e 857333, que contêm a posição do 
personagem na tela. 

Embora essa posição tenha sido co- 
locada nos endereços, ela ainda perma- 
nece no par HL. Todos os Id são essen- 
cialmente operações de cópia que colo- 
cam o valor na nova posição ou regis- 
tro, mantendo esse valor na posição ou 
registro de origem. No nosso caso, Isso 
é importante, pois, para chamar a roti- 
na de impressão de bloco em 58970, a 
posição de impressão na tela precisa es- 
tar no par de registros HL. 

O par BC, com o apontador de da- 
dos para a primeira figura de Willie sal- 
tando, é carregado em 57048. Esses da- 
dos foram apresentados em um artigo 
anterior de Avalanche e podem ser lidos 
numa tabela de dados na memória. 

A é carregado com 40 — código de 
azul sobre ciano, a cor de Willie — e 
DE, com 259. Isso faz com que a rotina 
de impressão de bloco imprima um blo- 
co um por três — 259 carrega | no re- 
gistro D e 3 no registro E; | x 256 + 3 
= 259. A rotina de impressão de bloco 
em 58970 é então chamada e imprime a 
figura de Willie saltando. 

Em seguida, o processador retorna da 
rotina. 
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WILLIE ESTÁ NO AR. 





A próxima parte da rotina de salto é 
chamada quando se aciona novamente 
a rotina de movimentação. Mais uma 
vez, essa seção começa verificando se a 
posição de memória 57335 contém 2 — 
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o que acontecerá se ela estiver sendo 
chamada pela segunda vez. 

Lembre-se de que o conteúdo de 
57335 está no acumulador quando o 
processador salta para essa parte da ro- 
tina; portanto, você não precisa carregá- 
lo de novo. Se o valor 2 está presente, 
o processador continua com a rotina e 
torna a incrementar o contador de sal- 
to com as instruções inc a e Id (57335),a. 
Assim, quando a rotina for chamada na 
próxima vez, a instrução cp 2 não irá en- 
contrar um valor 2 e a instrução jr 
nz,mjc mandará o processador para a 
terceira parte. 

Se o 2 estiver presente e esta parte da 
rotina for acionada, a nova posição de 
Willie no ar é carregada dos endereços 
57332 e 57333 para o par de registros 
HL. BC é carregado com 57000, que 
corresponde ao endereço inicial dos da- 
dos para a figura de Willie de pé com 
as pernas juntas. Esta é a figura que se- 
rá impressa quando, ao saltar, Willie es- 
tiver no alto. 

A é carregado com 40 — azul sobre 
ciano, a cor de Willie. DE é ajustado 
com 258 para imprimir um bloco de um 
por dois — | x 256 + 2 = 258. A roti- 
na de impressão em 58970 imprime Wil- 
lie no ar com as pernas juntas. 

O par DE é carregado com 64 e adi- 
cionado ao apontador de tela em HL — 
o que move esse apontador duas linhas 
para baixo. O par HL aponta agora pa- 
ra a cabeça de Willie. Lembre-se de que, 
subtraindo 64 (32 para cada linha), fa- 
zemos com que ele aponte para a posi- 
ção abaixo dos pés do personagem. A 
é carregado com 45 — ciano sobre cia- 
no, apenas um pedaço comum de céu. 
O apontador de dados BC é carregado 
com 15616, o início do conjunto de ca- 
racteres para um espaço em branco. 

A rotina print, no endereço 58217, é 
chamada, imprimindo um espaço em 


“branco para apagar a parte da figura de 


Willie que restou da primeira seção da 


| ed rotina de salto. Lembre-se de que um 





“bloco de um por três foi impresso ali, 


e, agora, apenas um bloco de um por 
dois. O processador retorna de novo, até 
a rotina voltar a ser chamada. 


DESCENDO 


Na próxima vez, o processador acio- 
nará a terceira parte da rotina de movi- 
mentação. Antes de mais nada, ele ve- 
rifica se é necessário acionar a quarta 
parte da rotina de salto. Em seguida, o 
contador de salto, que ainda está no 
acumulador, é incrementado e armaze- 


nado novamente em 57335, para que 0 


processador pule para a quarta parte da 
rotina na próxima vez. 

O apontador da posição de Willie na 
tela é carregado de volta em HL, O 
apontador de dados em BC é ajustado 
com o mesmo valor utilizado na primei- 
ra parte da rotina de salto. Willie des- 
cerá ao chão pelo mesmo processo que 
subiu. 

A é carregado de novo com 40 — 
azul sobre ciano; DE volta a ser ajusta- 
do com 259 — um por três. A rotina de 
impressão de bloco em 58970 é chama- 
da e imprime Willie descendo. 

O par DE é carregado com 32, valor 
adicionado ao par HL. O resultado é 
carregado de volta em 57332 e 57333. Is- 
so move o apontador de tela uma posi- 
ção para baixo. 

O processador retorna agora com as 
posições dos dados ajustados para a 
próxima parte. 


FIM DO SALTO 





Mais uma vez a rotina começa com 
um cp, para verificar se esta é a parte 
necessária da rotina de salto. Mas, ago- 
ra, estamos diante da última seção da 
rotina que faz Willie pular. Para que o 
processador chegue até aqui, o conteú- 
do do acumulador deve ser 4 (parte que 
finaliza o salto) ou 129, valor que indi- 
ca que Willie irá saltar para a frente. 
Nesse caso, jr nz,mfj manda o proces- 
sador para outra instrução de retorno, 
no final deste programa. O retorno se- 
rá apagado mais tarde, servindo para in- 
dicar o local onde se inicia a rotina de 
salto à frente, que será fornecida no pró- 
ximo artigo de Avalanche. Desta vez, 
entretanto, nenhum teste será feito e não 
precisaremos incrementar o contador de 
salto. É necessário que se ajuste o con- 
tador com 0, para que, quando a rotina 
de movimentação for novamente cha- 
mada, o processador execute a parte da- 
da no artigo anterior e verifique se al- 
guma tecla foi pressionada. Isso é feito 
por Id a,0 e ld (57335),a. 





A posição de Willie volta a ser carre- 
gada de 57332 e 57333 para o par de re- 
gistros HL, onde 32 é subtraído através 
do par DE, fazendo o apontador se mo- 
ver uma posição acima na tela. 

BC é carregado com o endereço de 
um espaço vazio e À, com o código de 
ciano sobre ciano. A rotina print, em 
58217, é chamada de novo. Apaga-se, 
assim, o caractere adicional ocupado 
por Willie, impresso na terceira parte da 
rotina de salto. 

A instrução jp 59153 faz o processa- 
dor voltar para o início da rotina de mo- 
vimentação, dada no artigo anterior. 
Como o contador de salto contém ago- 
ra o valor O, o processador executa a ro- 
tina que imprime Willie parado, apagan- 
do o que restou da antiga figura. 

A última instrução ret da listagem 
tem a função de prevenir que ocorra um 
erro no programa, caso a rotina de sal- 
to para a frente venha a ser chamada. 
Essa instrução será apagada na próxi- 
ma parte do programa. 


ERGUER SIE 


10 ORG 20140 

20 JUM JSR CLICK 
30 LDA 18261 

40 CMPA 41 

50 BNE MJA 

60 INC 18261 

70 LDX 18249 


BO PSHS X 

90 LEAX 256,X 
100 LDU 41536 
110 JSR CHARPR 
120 PULS X 

130 LEAX -256,X 
140 STX 18249 
150 LDU 417814 
160 JSR CHARPR 
170 LEAX 254,X 
180 LDU 417846 
190 JSR CHARPR 
200 RBRTS 

210 MJA CMPA 42 
220  BNE MJB 

230 INC 18261 
240 LDX 18249 
250 LEAX -256,X 
260 LDU 417870 
20 JSER CHARPR 
280 LEAX 254,X 
290  JSR CHARPR 
300 LEAX 254,X 
310 JSR CHARPR 
320 RTS 

330 MJB CMPA 43 
340  BNE MJC 

350 INC 18261 
360 LDX 18249 
380 LDU 41536 
390 JSR CHARPR 
400 LEAX 254,X 
410 LDU 417926 





LARS SA E E DNA 2 


420 JSR CHARPR 
430 LEAX 254,X 
440 LDU 417958 
450 JSR CHARPR 
460 LEAX 254,X 
470 LDU 417990 
480 JSR CHARPR 
490 ARTS 

500 MJC CMPA +4 
510 BNE MFJ 

520 CLR 18261 
530 LDX 18249 
540 PSHS X 

550 LDU 41536 
560 JSR CHARPR 
570 PULS X 

580 LEAX 256,X 
590 STX 18249 
600 LDU 417774 
610 JSR CHARPR 
620 LEAX 254,X 
630 JSR CHARPR 
640 ARTS 

650 MFJ RTS 

660 CHARPR EQU 19402 
670 CLICK EQU 20847 


A primeira coisa que a rotina JUM 
(ou de salto) faz é chamar a sub-rotina 
CLICK, que toca uma das notas do efei- 
to sonoro do salto. Se você ainda não 
tiver digitado essa rotina, adicione um 
retorno em 20847 por meio da instrução 
POKE, com 57 nesse endereço. Assim, 
não haverá erro quando o programa for 
testado. 





PULANDO POR ETAPAS... JÁ! 


O acumulador é carregado com o 
conteúdo da posição de memória 18261, 
que contém a variável do salto, Inicial- 
mente, ela está carregada com 1 ou 129 
— dependendo de qual tecla foi pressio- 
nada — indicando, portanto, se Willie 
deve pular verticalmente ou dar um sal- 
to à frente. Cada um desses tipos de pu- 
lo é dividido em quatro etapas, execu- 
tadas uma por vez, quando a rotina de 
movimentação é chamada. Para garan- 
tir o acionamento da etapa correta, a va- 
riável do salto é incrementada na exe- 
cução de cada parte. 

Se a variável do salto é O, utiliza-se 
a seção da rotina que examina o tecla- 
do para ver se uma ordem de pular foi 
dada. Caso a variável de movimentação 
contenha qualquer outro valor que não 
seja O, o processador vem para esta par- 
te do programa e executa uma das roti- 
nas de salto. 

Se for 1, a primeira parte da rotina 
é acionada. Durante sua execução, a va- 
riável de salto é incrementada, para que, 
na próxima vez que a rotina de movi- 
mentação for chamada, essa seção seja 
evitada e o processador passe para a par- 


te dois. Assim, sempre que a rotina de 
movimentação for acionada, ela chama- 
rá a próxima parte — até que a última 
(a quarta) seja executada e a variável do 
salto volte a ser ajustada com O. 

Se a última rotina ajustou a variável 
do salto com 129, esta parte é ignora- 
da. A rotina do salto para a frente, que 
daremos no próximo artigo de Avalan- 
che, é então acessada. 





O conteúdo do acumulador é compa- 
rado com 1. Se 1 não está presente, O 
processador passa para MJA. Caso con- 
trário, o processador continua com es- 
ta parte da rotina. 

A primeira coisa que ela faz é incre- 
mentar a variável em 18261, deixando- 
a pronta para a próxima vez. 

X é carregado com a posição de Wil- 


lie, que está no endereço 18249 e é ar- 


mazenado temporariamente na pilha. O 
apontador em X é incrementado com 
256, movendo-se um caractere para bai- 
xo. U é carregado com 1536, o endere- 
ço do canto superior esquerdo da tela. 
O processador chama então a rotina 
CHARPR, que imprime um bloco de 
céu azul sobre a metade inferior de Wil- 
lie. Se não fizermos isso, as pernas de 
Willie aparecerão em sua posição an- 
terior. 

A posição de Willie é novamente 
guardada na pilha e movida um carac- 
tere para cima pela subtração de 256. O 
resultado é armazenado de volta em 
18249, U é carregado com 17814, ende- 
reço inicial para os dados de Willie com 
as pernas abertas. Em seguida, a rotina 
CHARPR é chamada para imprimir a 
metade superior. 

O apontador de tela é incrementado 
com 254, passando a apontar para o iní- 
cio da metade inferior de Willie. U é re- 
carregado e a rotina CHARPR volta a 
ser chamada para imprimir sua metade 
inferior. 

O processador retorna depois de ter 
imprimido Willie com as pernas abertas 
um caractere acima do chão. 


Quando a rotina de movimentação de 
Willie for novamente chamada, o núme- 
ro no endereço 18261 será 2. Portanto, 
o processador irá pular a primeira par- 
te da rotina publicada no artigo anterior 
e a parte desta rotina que foi dada até 


agora. Mas, quando ele saltar de novo 
para MJA, o conteúdo de 18261, que 


ainda está no acumulador, será compa- 
rado com 2. Na chamada seguinte da ro- 
tina de movimentação, quando o con- 
teúdo do endereço 18261 tiver sido in- 
crementado mais uma vez, a instrução 
BNE MJB mandará o processador ve- 
rificar se seu valor é 3, 4, 129, 130, 131 
ou 132. Desta vez ele executará a próxi- 
ma seção da rotina, imprimindo a figu- 
ra de Willie no segundo estágio do pulo. 

Essa seção também começa incre- 
mentando o conteúdo de 18261, acer- 
tando-o para a próxima chamada da ro- 
tina de salto. Depois, a nova posição de 
Willie é carregada de 18249 para X e de- 
crementada de 256 — o que faz o per- 
sonagem se mover um caractere acima 
na tela. Mas agora o resultado não é co- 
locado de volta em 18261. Willie está no 
alto de seu pulo. 

U é carregado com 17870, endereço 
inicial dos dados da nova figura de Wil- 
lie. Desta vez, ele tem três caracteres de 
altura, impressos em três posições — 
CHARPR precisa ser chamada três ve- 
zes e X é decrementado com 254 para 
mover o apontador uma linha abaixo 
em cada chamada. 

Feito isso, o processador retorna. 





Quando a rotina de movimentação 
volta a ser chamada, o conteúdo de 
18261 é€3 e o processador executa a ro- 
tina MJB. Ela começa, mais uma vez, 
incrementando 18261, carregando 18249 
em X e subtraindo 256. 

Só que agora a figura anterior de 
Willie deve ser apagada. U é carregado 
com 1536, o endereço do canto superior 
esquerdo da tela, e a rotina CHARPR 
é chamada, imprimindo a cabeça de 
Willie na tela. O apontador X é incre- 
mentado para indicar a posição inferior 
na tela e U é carregado com o endereço 
dos dados para a impressão do resto da 
figura. 





QUATRO! 


Se nosso personagem está pulando na 
vertical, mais cedo ou mais tarde o pro- 
cessador irá executar MJC. As instru- 
ções CMPA 44 e BNE chamarão MJF 
apenas quando Willie for saltar na dia- 
gonal. Mas aquela pequena rotina será 
sempre executada. 

Como chegamos à quarta e última 
parte da rotina que faz Willie pular ver- 
ticalmente, o conteúdo da posição 18261 
não é incrementado. A instrução CLR 
18261 ajusta essa variável com O; a ro- 
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tina de movimentação de Willie é, por- 
tanto, executada e verifica-se, pelo exa- 
me do teclado, se ele deve pular de novo. 

A posição de Willie é carregada em 
X e colocada na pilha de máquina para 
armazenamento temporário. Como a fi- 
gura tem apenas dois caracteres de al- 
tura quando está parada ou andando, 
e três, quando está pulando, o caracte- 
re mais alto tem que ser apagado de no- 
vo. Note que, agora, Willie está numa 


posição intermediária, usada para tor- 
nar o pulo mais suave. 

U é carregado com 1536, que aponta 
para uma parte do céu no topo da tela. 
Quando é chamada, a rotina CHARPR 
apaga a cabeça da figura anterior. 

A posição de Willie é recuperada da 
pilha e adicionada a 256, o que faz a fi- 
gura se mover um caractere abaixo na 
tela. Lembre-se de que sua posição ti- 
nha sido decrementada com 256 no iní- 









cio desta rotina, quando Willie começou 
o salto. Assim, quando o novo valor é 
armazenado em 18261, o personagem 
simplesmente volta para a posição an- 
terior ao do pulo. 

U é então carregado com 17774, en- 
dereço inicial dos dados para a figura de 
Willie com as pernas juntas. CHARPR 
é chamada duas vezes e X é adicionado 
com 254, imprimindo a figura comple- 
ta de Willie na tela. 





"4 DEE RSENSE SS 


10 org 54854 

20 pi ld a, (-5202) 
30 Cp À 

40 jr nz,pb 

50 inc a 


60 ld (-5202),a 
70 ld h1,(-5205) 
BO ld de,32 
90 sbc hl,de 


100 ld (-5205),hl 
110 ld de, (62407) 
120 add hl,de 
130 ld a,1l0 

140 push hl 

150 call 77 

160 pop hl 

170 ld de,32 

180 add hl,de 
190 ld a,ll 

200 push hl 

210 call 77 

220 pop hl 

230 ld de,32 

240 add hl,de 
250 ld a,255 

260 call 77 

270 ret 

280 pb cp 2 

290 yr nz,pc 

300 inc a 

310 ld (-5202),a 
320 ld hº,(-5205) 
330 ld de, (62407) 
340 add hl,de 
350 ld a,0 

360 push hl 

370 call 77 

380 pop hl 

390 ld de,32 

400 add hl,de 


RN 5 a cm a AEE 


410 ld a,l 

420 call 77 

430 ret 

440 pc cp 3 

450 jr nz,pd 

460 inc à 

470 ld (-5202),a 
480 ld h1l,(-5205) 
490 ld de, (62407) 
500 add hl,de 
510 ld a,l0 

520 push hl 

530 call 77 

540 pop hl 

550 ld de,32 

560 add hl,de 
570 ld a,ll 

S80 call 77 

590 ret 

600 pd cp 4 

610 jr nz,mp 

620 ld a,0 

630 ld (-5202),a 
640 ld hi, (-5205) 
650 push hl 

660 ld de, (62407) 
670 add hl,de 
680 ld a,255 

690 call 77 

700 pop hl 

710 ld de,32 

720 add hl,de 

730 ld (-5205),h1 
740 jp 54603 

750 mp ret 

760 end 


A rotina de salto é dividida em qua- 
tro partes; cada uma delas é executada 
separadamente quando a rotina de mo- 
























vimentação de Willie é chamada. No úl- 
timo artigo, vimos como o processador 
é direcionado para a rotina pl quando 
o conteúdo da posição de memória 
— 5202 não é zero. Vimos também co- 
mo examinar o teclado para descobrir 
se as teclas de controle M e N foram 
pressionadas, individualmente ou ao 
mesmo tempo. 

Se a tecla M foi pressionada e Willie 
deve saltar, a posição de memória 
— 5202 é carregada com 1. Isso fará com 
que o processador seja direcionado pa- 
ra esta rotina na próxima vez que a ro- 
tina principal do jogo chamar a rotina 
de movimentação. 

Se Me N foram pressionadas, Willie 
deve saltar na diagonal e 129 é carrega- 
do em — 5202. Este caso será examina- 
do no próximo artigo. 

Por enquanto, suponhamos que a ro- 
tina de movimentação detectou que a te- 
cla N foi pressionada na última vez que 
foi chamada. Agora, ela voltou a ser 
chamada, o que fez o processador dire- 
cionar-se para a rotina do salto, cuja 
execução se inicia... 

O conteúdo da posição -— 5202 da 
memória é carregado novamente no acu- 
mulador. Em seguida, o valor em A é 
comparado com 1. Como a suposição 
feita exige que o valor seja 1, a instrução 
seguinte, jr nz,pb, não é executada e O 
processador continua com inc a. Com is- 
so, o conteúdo do acumulador passa a 
ser 2 — valor que é colocado de volta 
em — 5202 pela instrução ld ( — 5202), a. 
Incrementa-se, assim, o contador de sal- 
to para que, na próxima vez que a rotina 
de movimentação for chamada, ela acio- 
ne a rotina pl. Neste ponto, porém, a pri- 
meira parte do salto é deixada de lado 
e o processador executa a segunda. 

A posição de Willie na tela é carre- 
gada dos endereços — 5205 e — 5204 no 
par HL. O número 32 é subtraído dessa 
posição e colocado em DE; a instrução 
sbec h1, de é executada. Isso coloca o 
apontador de posição um caractere aci- 
ma de Willie na tela. O resultado é car- 
regado de volta nos endereços que 

contêm a posição de Willie, — 5205 

e — 5204. 

Embora esse valor tenha sido 
colocado nos endereços, ele per- 
manece no par HL. Todos os ld 

são essencialmente operações de 
cópia, que colocam o valor na nova 
posição ou registro, mantendo esse va- 
lor na posição ou registro de origem. No 
nosso caso isso é importante, pois o va- 
lor em HL é somado ao endereço inicial 
da Tabela de Nomes da VRAM, passan- 
do a apontar para a posição de Willie 
na TN. 

O acumulador é carregado com 10, 





IR 


código do primeiro padrão da figura de 
Willie subindo, com as pernas abertas. 
O apontador no par HL é guardado na 
pilha e a rotina 77 da ROM é chamada. 
Essa rotina escreve o valor de A no en- 
dereço da VRAM apontado por HL — 
o que equivale, aqui, a colocar o padrão 
na tela. 

O apontador é recuperado da pilha 
e somado com 32, voltando-se para uma 
posição abaixo da anterior. A é carre- 
gado com 11, código do segundo padrão 
da figura de Willie. O apontador é no- 
vamente guardado na pilha e a rotina 77 
da ROM é chamada. 

A figura de Willie saltando já está 
completa. Falta só apagar uma parte da 
figura anterior. Para isso, recupera-se O 
apontador da pilha e adiciona-se 32, fa- 
zendo-o apontar para o padrão que se 
pretende eliminar. O código do padrão 
de céu, 255, é colocado em A e a rotina 
77 é chamada, apagando parte da figu- 





ra anterior. Em seguida, o processador, 


retorna. 


“WILLIE ESTÁ NO AR 





A segunda parte da rotina de salto é 
chamada assim que se aciona a rotina 
de movimentação. Como as outras par- 
tes, ela começa verificando se a posição 
de memória — 5202 contém 2 — o que 
acontece se esta parte ainda não foi exe- 
cutada. 

Lembre-se de que o conteúdo de 
— 5202 está no acumulador quando o 
processador aciona esta parte — portan- 
to, não é necessário carregá-lo outra vez. 
Se o valor 2 está presente, o processa- 
dor continua com a rotina e incrementa 
o contador de salto com as instruções 
inc a e ld (— 5202), a. Quando esta par- 
te for acessada na próxima vez, cp 2 não 
irá encontrar o valor 2 e jr nz,pc man- 
dará então o processador para a ter- 
ceira parte do salto. 

Se o valor 2 estiver presente e esta 
parte da rotina foi acionada, a atual po- 
sição de Willie é transferida dos ende- 
reços — 5205 e — 5204 para o par de re- 
gistros HL. O par DE é carregado com 
o endereço inicial da TN da VRAM e so- 
mado em HL. O acumulador é carrega- 
do com 0, o código do primeiro padrão 
da figura de Willie parado. O aponta- 
dor em HL é guardado na pilha e a ro- 
tina 77 da ROM é chamada. 

O apontador é recuperado da pilha 
e somado com 32 por meio de DE, o que 
o faz se mover uma posição para bai- 
xo. A é carregado com 1, o código do 
segundo padrão da figura. A rotina 77 
é chamada de novo. Neste ponto, Wil- 
lie está impresso no ar com as pernas 
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juntas. Observe que agora não é preci- 
so apagar nada, pois a figura foi impres- 
sa sobre a anterior. 

O processador retorna da rotina até 
a próxima chamada. 


DESCENDO 





Desta vez, O processador executa a 
terceira parte da rotina de salto. Primei- 
ro, ele verifica se é necessário acionar a 
quarta parte do salto. 

Em seguida, o contador de salto, que 
ainda está no acumulador, é incremen- 
tado e devolvido a — 5202, para que o 
processador passe para a quarta parte 
na próxima vez. 

O apontador de Willie na tela é car- 
regado no par HL. O par DE é carrega- 
do com o endereço inicial da TN da 
VRAM e somado em HL. Em seguida, 
a mesma figura da primeira parte é im- 
pressa de maneira idêntica, utilizando a 
rotina 77 da memória ROM. O proces- 
sador retorna desta parte. 





FIM DO SALTO 


Mais uma vez a rotina inicia com um 
cp, verificando se esta é a parte neces- 
sária da rotina de salto. Para que o pro- 
cessador chegue até aqui, o conteúdo do 
acumulador deve ser 4 ou 129. Se for 4, 
esta é a última parte da rotina do salto; 
se for 129, Willie irá saltar na diagonal. 
Nesse caso, a instrução jr nz,mp man- 


da o processador para a instrução de re- 
torno do final deste programa. Essa ins- 
trução será apagada pela rotina de sal- 
to diagonal, que apresentaremos no pró- 
ximo artigo de Avalanche. 

Nenhum teste será feito, nem preci- 
saremos incrementar o contador de sal- 
to. É necessário, porém, ajustar esse 
contador com zero, para que, quando ' 
a rotina de movimentação for chamada 
de novo, o processador execute a primei- 
ra parte, dada no artigo anterior, que 
verifica se alguma tecla foi pressionada. 
Isso é feito pelas instruções Id a,0 e ld 
(— 5202),a. 

A posição de Willie é mais uma vez 
carregada no par HL, sendo armazena- 
da na pilha. O endereço inicial da TN 
da VRAM é colocado em DE e somado 
ao valor de HL. O acumulador é carre- 
gado com 255, o código do padrão de 
céu. A rotina 77 da ROM é chamada, 
apagando parte da figura anterior. 

A posição de Willie é recuperada da 
pilha para HL, somada com 32 através 
de DE e, finalmente, devolvida aos en- 
dereços — 5205 e — 5204. 

A primeira rotina de movimentação 
de Willie é chamada por jp 54603. Não 
se preocupe com o que restou de sua fi- 
gura, pois nossa rotina apagará todos os 
vestígios indesejáveis. 

O último ret da listagem é colocado 
apenas para prevenir possíveis erros, que 
ocorreriam se a rotina de salto diagonal 
fosse chamada sem estar presente na me- 
mória. Essa instrução será apagada pe- 
la parte do programa dada no próximo 


“* artigo de Avalanche. 
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MODELOS 





DA REALIDADE 


Os engenheiros constroem modelos 
em miniatura, para testar suas proprie- 
dades. Quando fazem isso, estão simu- 
lando um sistema. O programa do arti- 
go da página 1121 executa algo bem pa- 
recido: ele simula um evento, oferecen- 
do resultados compatíveis com a reali- 
dade. Veremos agora como combinar os 
princípios da simulação e da estatística 
para resolver problemas de ordem prá- 
tica e para acrescentar a nossos jogos 
uma dose de realidade. 


| TIPOS DE MODELO 


Um modelo é uma representação sim- 
plificada de um sistema real, por meio 
da qual se procura explicar ou prever as 
características e o comportamento des- 
se sistema. Existem três tipos de mode- 
lo — icônico, análogo e simbólico —, 
nem todos úteis aos usuários de micro- 
computadores. 

Modelos icônicos não possuem par- 
tes móveis. Maquetes de estradas € pré- 
dios incluem-se nessa categoria. Ao pro- 
jetar a construção de uma ponte sobre 
um estuário, por exemplo, engenheiros 
civis providenciam uma maquete da 
ponte e da área adjacente, a fim de es- 
tudar os efeitos da maré, da correnteza 
e do vento. Embora se utilizem compu- 
tadores para interpretar os resultados, 
a necessidade de coletar dados através 
de um experimento torna a simulação 
dispensável. 

Modelos análogos representam um 
fenômeno — ou quantidade — por ou- 
tro. É o caso de um economista que, pa- 
ra representar o fluxo de dinheiro na 
economia mundial, utiliza o fluxo de 
água de um conjunto de tanques cheios 
até diferentes níveis. Esse tipo de repre- 
sentação pode ser feito eletronicamente 
— empregando-se computadores analó- 
gicos e não digitais. 

O terceiro tipo de modelo, o simbó- 
lico, é o que apresenta particular inte- 
resse para o usuário do micro. 








INCERTEZAS 


Num modelo simbólico, representa- 
se determinado sistema por meio de uma 
























fórmula matemática. Para estudar a dis- 
tância S percorrida por um carro que se 
desloca à velocidade constante Y por um 
certo período de tempo T, um modelo 
adequado seria S = V*T. Essa equação 
é chamada determinística, pois não 
comporta nenhum elemento de incerte- 
Za OU acaso. 

Existem, contudo, eventos que não 
podem ser previstos com tal grau de pre- 
cisão. Temos um bom exemplo disso no 
programa apresentado no artigo da pá- 
gina 776, que simula o lançamento de 
uma moeda. Modelos que descrevem fe- 
nômenos que incluem um elemento de 
acaso são chamados estocásticos. Este 
é o tipo de evento que se costuma simu- 
lar, especialmente em jogos. 

Algumas variáveis só podem assumir 
valores discretos como 0, 1, 2, 3... O nú- 
mero de gols em uma partida de fute- 
bol pertence a essa categoria. Quando 


à “—=— as, = E ai = - si 5 as is od 
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Para ganhar a quina da loto, basta 
escolher cinco números ao acaso — 
os cinco números certos, é claro. 
Veja aqui algumas dicas estatísticas 
para programar eventos aleatórios. 


simulamos variáveis aleatórias como es- 
sa, precisamos sempre levar em conta 
que certos processos aparentemente di- 
ferentes têm, na realidade, a mesma es- 
trutura. 

Tomemos como exemplo o programa 
que simula o lançamento da moeda. Ele 
poderia ser facilmente modificado a fim 
de fornecer os resultados de um jogo de 
dados ou da cobrança de pênaltis. Em 
todas essas circunstâncias há um certo 
número de eventos independentes, cada 
qual com uma probabilidade definida. 
No caso de lançarmos uma moeda qua- 
tro vezes, vários resultados são possi- 
veis: quatro caras, três caras e uma co- 
roa, duas caras e duas coroas, uma ca- 
ra e três coroas, quatro coroas. Cada um 
dos resultados, porém, tem uma proba- 
bilidade definida. Os processos desse ti- 
po, conhecidos como processos de Ber- 
noulli, são formados por variáveis alea- 
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tórias discretas. Eles não se prestam, por 
exemplo, para descrever o número de 
acidentes do campeonato mundial de 
Fórmula-l ou o número de chamadas te- 
lefônicas que chegam a uma central em 
meia hora. Nessas situações, OS eventos 
ocorrem ao acaso durante um certo pe- 
ríodo de tempo e não são o resultado de 
experimentos sucessivos. Tais processos 
são denominados processos de Poisson. 

Para obter uma distribuição de Pois- 
son, digite este programa: 


10 BORDER O: 
CLS 

20 POKE 23658,0 
40 PRINT AT 1,6; 
IMULACAO DE POISSON *: 
: PRINT 

50 INPUT "QUAL E O VALOR MEDI 


PAPER O: INK 7: 


INVERSE 1;"5 
PRINT 





1 TIPOS DE MODELO 
[E INCERTEZAS 
a PROBABILIDADE E 

COMPORTAMENTO ALEATÓRIO 
[e 





VARIÁVEIS ALEATÓRIAS 





O ":a 
70 INPUT "TAMANHO DA AMOSTRA 


BO FOR 1=1 TO n 

90 LET c=0: LET t=1 

LET s=EXP (-a) 

LET t=t* (RND*1) 

IF t<=s THEN PRINT ” 
"scr: GOTO 150 

130 LET c=c+1 

140 GOTO 110 

150 NEXT à 

160 INPUT ” OUTRA AMOSTRA 

N) ?7":95 

170 IF qg58="s" THEN GOTO 10 
180 STOP 


io 


10 R=BND(-TIME) 

30 CLS:WIDTH 40 

40 PRINT "Distribuão de Poisson 
"PRINT: PRINT 

50 INPUT "Qual é o valor da méd 
1a" ;Ã 

60 PRINT 

70 INPUT "Tamanho da amostra”;N 
BO FOR I=1 TON 

90 C=0:T=1 

100 S=EXP(-A) 

110 T=T*RND(L) 

120 IF T<=S THEN PRINT LEFTS(ST 
ARS(C)+" ",8); :GOTO 150 
130 C=C+1 

140 GOTO 110 
150 NEXT I 

160 PRINT: INPUT 
(Y/N)";6S 

170 IF G8="Y” 
180 END 


30 HOME 

40 PRINT TAB( 9);"DISTRIBUICA 
O DE POISSON": PRINT PRINT 

S0 INPUT "QUAL O VALOR DA MEDI 
A "ih 


(5/ 


"Outra amostra 


THEN 30 


60 PRINT 

70 INPUT "TAMANHO DA AMOSTRA ” 

:N 

B0 FOR I =1TON 

90 C=0:T=1 

100 S = EXP ( - A) 

110 T=TÃx* BND (1) 

120 IF T< = S THEN PAINT XL 

EFTS ( STARS (C) + * "28)1 
Goto 150 

130 C=C+1 

140 GOTO 110 

150 NEXT I 

160 PRINT INPUT "OUTRA AMOST 

RA (S/N) ";G58 





| “SIMULAÇÕES 
o DISTRIBUIÇÃO UNIFORME 
B DISTRIBUIÇÃO NORMAL 
EE DISTRIBUIÇÃO EXPONENCIAL 
[E COMPARAÇÃO DE EVENTOS 


170 THEN 30 


180 


30 CLS 

40 PRINT €6,"simulacao de poiss 
on” : PRINT: PRINT 

50 INPUT"QUAL E O VALOR MEDIO ” 
:A 

60 PRINT 

70 INPUT"TAMANHO DA AMOSTRA 
BO FOR I=1l TON 

90 C=0:T=1 

100 S=EXP(-A) 

110 T=T*RND(O0O) 

120 IF T<=S THEN PRINT LEFTS(ST 
AS (C)+" ",B);:GOTO 150 
130 C=C+1 

140 GoToO 110 

150 NEXT 1 

160 PRINT: INPUT"OUTRA AMOSTRA ( 
S/N) ";G5 

170 IEF G8="S” THEN 30 

180 END ; 


IF 68 = 5” 
END 


*+*N 


Ao ser executado, o programa soli- 
cita o valor médio, bem como o núme- 
ro de elementos da amostra. A parte 
principal do programa (linhas 80 a 150) 
Usa uma equação para gerar as variáveis 
de Poisson. A variável S recebe o valor 
de e (uma constante matemática) eleva- 
da ao valor de A (valor da média). A li- 
nha 110 escolhe um número aleatório 
entre 0 e T. A linha 120 compara as va- 
riáveis Se T para decidir se imprime ou 
não uma variável €. 

Suponhamos que você queira criar 
um jogo onde uma nave espacial enfren- 
ta uma chuva de meteoritos. Se o núme- 
ro médio de “colisões” em um período 
de um minuto é dois, e você precisa de 
cinco simulações de um minuto, o com- 
putador poderá fornecer o resultado 2, 
3, 2, 1, O como amostra do número de 
choques que a nave sofrerá em um mi- 
nuto. Cada valor não é muito diferente 
de 2 (a média) e há cinco valores — o 
tamanho da amostra. 

Poderíamos usar a mesma técnica pa- 
ra simular os resultados de uma série de 
partidas de futebol? É claro que sim. Es- 
taremos lidando com uma variável dis- 
creta, pois o número de gols em uma 
partida só pode assumir valores inteiros. 
Gols são eventos que ocorrem ao aca- 
so, em determinado espaço de tempo — 
e não o resultado de um certo número 
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de experimentos repetidos. Portanto, 
nesse caso, um processo de Poisson é 
adequado. 

Além do número de partidas, preci- 
saremos escolher a média de gols por 
partida. Este é o momento de buscar da- 
dos no mundo real. A tabela abaixo 
mostra a média de gols por partida nos 
campeonatos da primeira divisão ingle- 
sa entre 1977 e 1983. 





Baseados nesses valores, podemos es- 
timular uma média de 1.7 gols por par- 
tida para o time da casa e de 1.0 para 
o time visitante. Dividindo esses núme- 
ros por dois, obteremos as médias de 
0.85 e 0.5 para meio tempo de jogo. De 
posse desses dados, fica fácil elaborar 
o programa: 


20 DIM ag(25,10): DIM h5(25, 
10): DIM h(50): DIM a(50): 
DIM £(50): DIM q(50) 


30 BORDER O: PAPER O: INK 7: 
CLS 
50 PRINT AT 0,8; INVERSE 1;” 


PLACAR FINAL *” 

70 INPUT * QUANTAS PARTIDAS ( 
I=22) ? “qn 

75 IF n<l OR n>25 THEN GOTO 
70 
B0 FOR i=1 TO n 
90 PRINT "PARTIDA" ;1 


100 INPUT "TIME DA CASA”;hS(i) 
110 INPUT "TIME VISITANTE” ;as( 
1) 

120 NEXT à 

130 PAUSE 400 

140 PRINT 

150 FOR i=1l TO 2*n 

160 LET c=0: LET t=1 

170 LET s=EXP (-.85) 

180 LET t=t*(RND*1) 

190 IF t<=s THEN LET h(i)j=c: 
GOoTo 220 

200 LET c=c+1 

210 GOTO 180 

220 NEXT i 

230 FOR i1=1 TO 2*n 

240 LET c=0: LET t=1 

250 LET a=EXP (-.5) 


260 LET t=t*(AND*1) 
270 IF t<=s THEN LET a(i)=c: 
GoTo 300 | 


280 
290 
300 


LET c=ctl 

GoTo 260 

NEXT à 

310 CLS 

320 PRINT TAB (7);" PLACAR PRIM 
EIRO TEMPO" ;'** 

330 FOR i=1 TO n 

340 PRINT hS(i);h(i); TAB 20;aS 
(1);a(i): PRINT 

350 NEXT à 

360 PRINT ” QUALQUER TECLA PAR 
à CONTINUAR” 

370 IF INKEY5S="" THEN GOTO 


370 

375 CLS 

380 PRINT TAB 12;"PLACAR FINAL 
390 FOR i=1 TO mn 
400 LET £f(i)=h(i)+h(n+1) 
410 LET glij=a(i)ta(i+n) 
420 PRINT hS(i);f(i);TARB 20;a5 
(i);o(i): PRINT 
430 NEXT à 
440 INPUT ” NOVAMENTE ? (m/nm 
) “198 
450 IF g8="n”" THEN GOTO 490 
460 INPUT "MESMOS TIMES (s/n) 
?" :p$8 


470 IF p$="n”" THEN GOTO 70 
480 GOTO 150 
490 STOP 


5 BR=ARND(-TIME) 

20 DIM AS(25),H5(25),H(50),A(50 
3), FA(25),FH(25) 

30 CLS:WIDTH 40 

50 PRINT TAB(13);"SCORE FINAL” 
60 PRINT:PRINT 


70 INPUT"Quantas partidas (1-25 
)74N 

75 IF N<l OR N>25 THEN 70 

BO FOR I=1l TON 

90 PRINT "PALITO";I 

100 INPUT "Time da casa" ;HS(T) 
110 INPUT "Time de fora";AS(TI) 
120 NEXT IT 

130 FOR V=1 TO 2000:NEXT V 

140 PRINT 


150 FOR I=1 
160 C=0:T=1 
170 S=EXP(-.85) 

180 T=T*RND (O) 

190 IF T<=S8 THEN H(I)=C:GOTO 22 
O 

200 €C=C+1 

210 GOTO 180 

220 NEXT I 

230 FOR I=1l TO 2*N 

240 C=0:T=1 

250 S5=EXP(-.5) 

260 T=T*RND(0O) 

270 IF T<=S8 THEN A(I)=C:GOTO 30 
O 

280 C=C+1 

290 GOTO 260 

300 NEXT I 

310 CLS 

320 PRINT TAB(9);"Score temporá 
rio”: PRINT: PRINT 

330 FOR I=1 TON 

340 PRINT HS(I):;H(I) ASCIJGA(I) 
350 NEXT I 


TO 2*N 


360 
a” 
370 
380 
390 
400 
410 
420 
I) 
430 
440 
"4163 
450 
460 
PS 
a7o 


P..INT "Aperte qualquer fecl 


IF INKEYS="" THEN 370 


CLS 


FOR I=1 TON 
FH(I)=H(I)+H(N+I) 
FA(I)=A(I)J+A(N+I) 

PRINT HS(I);FH(I),AS(CI);FA( 


NEXT 
PRINT: INPUT"OUTRA VEZ (5/N) 


IF GS="N”" THEN END 
INPUT "Mesmos times (S/N)”"; 


IF P$="N” THEN 70 


480 GOTO 150 


Dat 


20 DIM AS(25),H$(25),H(50) ,A(S 
0) FA(25) , FH(25) 

30 HOME 

50 PRINT TAB( 13);"SCORE FINA 
| 

60 PRINT : PRINT 

70 INPUT "QUANTAS PARTIDAS (1- 
25) ";N | 
75 IFN<1O0RN>gZ5 THENJZO 

80 FORI=1TON 

90 PRINT "PARTIDA ";1 

100 INPUT "TIME DA CASA ";HS(I 
) 

110 INPUT "TIME DE FORA ";AS(I 
) 

120 NEXT 1 

130 FOR V = 1 TO 2000: NEXTV 
140 PRINT 

150 FORI =1TO2Z4%*N 

160 C=0O:T=1 

170 8 = EXP ( - .85) 

180 T = T * BND (1) 

190 IF TT< =STHENH(I)=c 
+ GOTO 220 

200 C=C+1 

210 GOTO 180 

220 NEXT 1 

230 FORI =1TO2Z*N 

240 C= 0:T=1 

250 S = EXP ( - .5) 

260 T = Tx* BRND (1) 

270 IF T< =STHENA(I) =c€: 
GOTO 300 

280 C=C+1 

290 GOTO 260 

300 NEXT I 

310 HOME 

320 PRINT "SCORE PARCIAL ": PR 
INT : PRINT 

330 FORI =1 TON 

340 PRINT HS(I);H(I),AS(I):A(I 
) 

350 NEXT 1 

360 PRINT "APERTE QUALQUER TEC 
LA * 

370 GET G8 

375 HOME 

380 PRINT TAB( 13);"SCORE FIN 
AL”: PRINT : PRINT 

390 FORI =1 TON 

400 FH(I) = H(I) + H(N + 1) 

410 FA(I) = A(I) + A(N + 1) 

420 PRINT H$S(I);FH(I),AS(I);FA 
(1) 
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430 NEXT 

440 PRINT : INPUT "OUTRA VEZ ( 
S/N) ":PS 

450 IF PS = "Nº THEN END 

460 INPUT "MESMOS TIMES (S/N) 
“168 

470 IF GS = "Nº THEN 70 

480 GOTO 150 


20 DIM AS(25),H8(25),H(50),A(5SO 
) FA(25), FH(25) 

30 CLS 

50 PRINT €11,"placar final” 

60 PRINT:PRINT 

70 INPUT"QUANTAS PARTIDAS (1-25 
) TN 

BO FOR I=1l TON 

90 PRINT"PARTIDA ";I 

100 INPUT"TIME DA CASA ";HS5(I) 

110 INPUT"VISITANTE ";AS(I):PRI 
NT 

120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
2/0 
280 
290 
300 


NEXT I 
FOR V=1 
PRINT 
FOR I=1 TO 
C=0:T=1 
S=EXP(-.B5) 
T=TX*AND (0) 
IF T<=S THEN H(1)=C:GOT0O220 
C=C+1 

GOTO 180 

NEXT I 

FOR 1I=1 TO 2*N 

C=0:T=1 

S=EXP(-.5) 

T=T*END (0) 

IF T<=S5 THEN A(I)=C:G0T0O300 
E=C+1 

GOTO 260 

NEXT 1 

310 CLS 

320 PRINT €65,"placar primeiro t 
empo” : PRINT:PRINT 

330 FOR I=1 TON 

340 PRINT HS(I);H(I)J;ASCI)JGA(I) 
350 NEXT 1 
360 PRINT” 
CONTINUAR” 
370 IF INKEY5S="" THEN 370 

375 CLS 

380 PRINT €10,"placar final”:PR 
INT: PRINT 

390 FOR 1I=1 TON 

400 FH(I)=H(I)+H(N+1) 

410 FA(I)J=A(I)J+A(N+I) 

420 PRINT HS(I);FH(I),AS(I)GEA( 
1) 

430 NEXT 

440 PRINT: INPUT"NOVAMENTE 
) "4568 

450 IF GS="Nº THEN END 
460 INPUT"MESMOS TIMES 
:PS 

470 IF PS="N” 
480 GOTO 150 


Da linha 20 a 120 o programa solici- 
ta as informações iniciais. Para evitar 
muita digitação, podemos substituir os 
nomes dos times por letras. As linhas 
150 a 300 usam o algoritmo de Poisson, 


TO 2000:NEXT V 


2Z2*N 


QUALQUER TECLA PARA 


(S/N 


(S/N) * 


THEN 70 


idêntico ao do programa anterior, para 
gerar os resultados parciais e finais dos 
jogos. As linhas 320 a 350 cuidam da 
saída dos resultados parciais e as linhas 
390 a 430, dos resultados finais. 

A simulação de Poisson funciona 
bem com variáveis discretas. Algumas 
vezes, porém, precisamos utilizar variá- 
veis fracionárias ou continuas. Em uma 
simulação de competição de salto em 
distância, por exemplo, podemos obter 
uma medida qualquer entre sete e nove 
metros — o comprimento do salto é 
uma variável aleatória, continua. 

Embora a função RND seja adequa- 
da para simulações de Poisson com nú- 
meros inteiros, será necessário modificá- 
la um pouco para utilizá-la com variá- 
veis contínuas. O quadro da página 1180 
compara as distribuições obtidas com os 
dois tipos de variáveis: variáveis discre- 
tas resultam em distribuições uniformes, 
enquanto variáveis continuas resultam 
em distribuições exponenciais e distri- 
buições normais. 

Digite o próximo programa para 
transformar os números aleatórios uni- 
formemente distribuídos, criados pela 
função RND do computador, em núme- 
ros distribuídos exponencialmente. 


10 BORDER O: INK 7: PAPER O: 
CLS 

20 POKE 23658,0 

30 PRINT AT 0,4; INVERSE 1;" 


SIMULACAO EXPONENCIAL *'** 
50 INPUT "QUAL E O VALOR MEDI 


O? Tia 
70 INPUT "TAMANHO DA AMOSTRA 
7” “um 
80 FOR i=1 TO nm 
90 LET x=(-a) *LN (RND*1): LET 
vy=INT (x*10): PRINT * a 
1*Y, 
100 NEXT à | 
110 INPUT ” OUTRA SIMULACAO (Bm 
in) ? "i;98 
120 IF q8="a” THEN GOTO 30 
130 STOP 


py 


10 R=RND(-TIME) 
30 CLS:WIDTH 40 
40 PRINT TAB(B);"Simulação expo 
nencial”:PRINT:PRINT 
50 INPUT "Qual o valor da média 
":h 
60 PRINT 
70 INPUT "Tamanho da amostra";N 
BO FOR 1£=1 TON 
90 X=(-A) *LOG (BND (1)) :T=INT(X*1 
0) :PRINT LEFTS (STR$(Y/10)+" 

" 8); 
100 NEXT I 
110 PRINT: INPUT "Outra vez (S/N 
)":68 


120 IF G3="8" THEN 30 
130 END 


Lob] [çó] 


30 HOME 

40 PRINT TAB( 7);"SIMULACAO E 
XPONENCIAL": PRINT : PRINT 

50 INPUT "QUAL O VALOR DA MEDI 
A "IA 





60 PRINT 

70 INPUT "TAMANHO DA AMOSTRA " 
:N 

BO FOR I = 1TON 

90 X = ( - A) * LOG ( BND (1)) 
To= INT (X * 10): PRINT LEFT 
S ( STR$S (Y / 10) + * "8 
Ji 

100 NEXT 1 

110 PRINT : INPUT "OUTRA SIMUL 
ACAO (8/N) ";G38 

120 IF G$ = "8" THEN 30 

130 END 

30 CLS 


40 PRINT €5,"simulacao exponenc 

ial”": PRINT:PRINT 

50 INPUT"QUAL E O VALOR MEDIO ” 

A 

60 PRINT 

70 INPUT"TAMANHO DA AMOSTRA ";N 

BO FOR I=1 TON 

90 Xx=(-A) *LOG(RND(0)) :Y=INT(X*1 

0) :PRINT LEFTS(STRS(Y/10)+" 
"2.8)5 

100 NEXT I 

110 PRINT: INPUT"OUTRA SIMULACAO 
(S/N) ";65 

120 IF G$="5" THEN 30 

130 END 


Devemos informar ao programa, as- 
sim que o executamos, O valor da mé- 
dia. Feito isso, os valores da amostra 
são exibidos. A linha 90 realiza toda a 
tarefa, utilizando um recurso matemá- 
tico chamado método da transformação 
inversa, para mudar a distribuição das 
variáveis. Esse programa pode ser usa- 
do para simular, por exemplo, o tempo 
de espera em um lava-rápido ou a du- 
ração das reservas de combustível em 
um jogo espacial. Em ambos os casos, 
o tempo decorrido é o elemento central. 





DISTRIBUIÇÃO NORMAL 


A curva normal é o tipo mais conhe- 
cido de distribuição estatística. Tem nu- 
merosas aplicações, descrevendo toda 
espécie de fenômeno natural — como a 
estatura e o peso de homens adultos. 
Além disso, constitui o meio mais apro- 
priado para a representação de uma dis- 
tribuição de erros. 

Como exemplo do uso da distribui- 
ção normal em um modelo, considere- | 
mos esta situação: as vendas (S) de uma | 





|| | 








empresa dependem de um valor-base 
(So) e do total gasto em propaganda (A). 
A relação poderia ser resumida na ex- 
pressão S = So + (B*A) + desvio nor- 
mal, onde B é uma constante conhecida. 

Digite o programa a seguir para ver 
como funciona. 


10 BORDER O: PAPER O: INK 7: 
CLS 

20 DIM u(50): DIM v(50): DIM 
z (50) 

50 PRINT AT 0,7; INVERSE 1; 


SIMULACAO NORMAL "*" 

60 INPUT ” QUAL E O VALOR MED 
IO 7? F:a 

70 INPUT ” QUAL E O VALOR MIN 
IMO ? F;m 

BO INPUT ” TAMANHO DA AMOSTRA 
? “uam 

90 LET b=(a-m)/2.5 


1=1 TO n 

t=0 

3=1 TO 15 

Y (3) =BND*1 

LET t=t+y (53): NEXT 3 

150 LET u(i)=((t/15)-.5)/SQR ( 
1/(12*15)) 


FOR 
LET 
FOR 
LET 


160 LET z(ij=a+(b*u(i)) 

170 LET p=INT (10*z(1)): PRINT 
a "1.1*p, 

180 NEXT 1 

190 INPUT ” OUTRA AMOSTRA *? nd 
gs 

200 IF qg8="B” THEN GOTO 50 
210 STOP 


| pi 
a | 
e] 
Fi “di 
F Em " 


10 R=BND(-TIME) 

30 CLS:WIDTH 40 

40 DIM U(50),Y(50),Z(50) 

50 PRINT "Simulação normal”:PRI 
NT: PRINT 

60 INPUT "Qual o valor da média 
TI:Ã 

70 INPUT "Qual o valor mínimo”; 


M 

80 INPUT "Tamanho da amostra"”;N 
90 B=(A-M)/2.5 

100 FOR I=1 TON 

110 T=0 

120 FOR J=1 TO 15 

130 Y(J)=RND(1) 

140 T=T+Y (JJ) :NEXT J 

150 U(I)=((T/15)-.5)/S0QR(1/(12* 
15)) 

160 Z(I)=A+(B*U(I)) 

170 P=INT(10*Z(1)):PRINT LEFTS( 
STARS (P/10)+" PB) 

180 NEXT I 

190 PRINT: INPUT"Outra amostra”; 
GS 

200 IF G$="S" THEN 50 

210 END 


30 HOME 
40 DIM U(50),7(50),Z(50) 
50 PRINT TAB( 9);"DISTRIBUICA 


O NORMAL”: PRINT PRINT 

60 INPUT "QUAL O VALOR DA MEDI 
A "IA 

70 INPUT "QUAL O VALOR MINIMO 
":M 

80 INPUT "TAMANHO DA AMOSTRA " 
“N 


90 B= (A-M)/ 2.5 

100 FOR I = 1 TON 

110 T = 0 

120 FORJ = 1 'TOlS 

130 Y(J) = BND (1) 

140 T = T+ Y(J): NEXT J 
150 U(I) = (CT / 15) - «3) 
OR (1 / (12 * 15)) 


f Ss 


160 Z(1) = A + (B * U(I)) 
170 P = INT (10 * Z(1I)): PRINT 
LEFTS ( STR$S (P / 10) + 7 
" 18) 
180 NEXT 1 
190 PRINT : INPUT "OUTRA  AMOS 
TRA (S/N) ";G8 
200 IF GS = "8" THEN 50 
210 END 


30 CLS 

40 DIM U(50),Y7(50),2Z(50) 

S0 PRINT €7,"distribuicao norma 
1”: PRINT: PRINT 

60 INPUT"QUAL E O VALOR MEDIO ” 


70 INPUT"QUAL E O VALOR MINIMO 
80 INPUT"TAMANHO DA AMOSTRA "iN 
90 B=(A-M)/2.5 


100 FOR I=1 TO N 

110 T=0 ; 

120 FOR J=1 TO 15 

130 Y(J)=RND(O) 

140 T=T+Y (J) NEXT Jd 

150 U(I)=((T/15)-.5)/SQR(1/(12* 
15)) 


160 Z(I)=A+(B*U(I)) 
170 P=INT(10*Z(1)):PRINT 
STR$(P/10)+" "BJ; 


LEFTS( 


180 NEXT 1 

190 PRINT: INPUT"OUTRA AMOSTRA " 
:65 

200 IF G$="5" THEN 50 

210 END 


Execute O programa e forneça o va- 
lor médio e o valor mínimo, bem como 
o tamanho da amostra. Tecnicamente, 
o intervalo de uma distribuição normal 
é infinito, o que deixa uma chance mui- 
to pequena de se obter um valor menor 
que o mínimo especificado. 

“As linhas 120 a 150 usam a função 
RND para criar quinze números entre 
O e 1, somando-os em seguida. A linha 
160 calcula a média e o fator de escala. 
Depois de ter sido escalonado, o resul- 
tado é impresso. 

Para maior clareza, poderá ser útil 
criar duas distribuições com a mesma 
média e aproximadamente o mesmo in- 
tervalo, mas com formatos diferentes. 
Execute o programa e forneça os valo- 
res 100 para a média, 50 para o mínimo 
e 40 para o tamanho da amostra. Se os 
valores obtidos fossem plotados, teria- 
mos um gráfico semelhante ao da figu- 
ra €C (quadro da página 1180). 

Agora, apague a linha 90 e as linhas 
120 a 160. Substitua a linha 110 pela se- 
guinte: 


110 Z(I)=M+RND(1)*2* (A-M) 


No TRS-Color, use RND (0) no lu- 
gar de RND (1). 

Com essas alterações, a distribuição 
normal é convertida em distribuição uni- 
forme. Execute o programa novamente 
e entre os valores especificados, compa- 
rando os resultados. Desta vez, OS va- 
lores levariam a um gráfico como o da 
figura A (página 1180). 

Observe que a simulação normal acu- 
mula maior número de resultados pró- 
ximos à média. 





LINHA 


Apple Il + 

Apple Il+ 
Apple Il + 
Apple Il + 
Apple Il + 
Apple ll + 
Apple + 
Apple ll+ 
Apple Il + 
Apple + 
Apple + 
Apple + 
Apple ll+ 
Apple ll + 
Apple ll+ 
Apple ll + 
Apple Il + 
Apple ll + 
Apple Il + 
Apple + 
Apple Hl+ 
Apple Il + 

Apple Il + 
Apple ll + 
Apple ll + 
Apple ll + 
Apple Il + 
Apple ll+ 
Apple lle 

Apple Ile 


Apple Ile 
MSX 


MSX 


FABRICANTE Po 


Appletronica 
CCE 

CPA 

CPA 

Digitus 
Dismac 
ENIAC 
Franklin 
Houston 
Magnex 
Maxitronica 
Maxitronica 
Maxitronica 
Maxitronica 
Microcraft 
Milmar 
Milmar 
Milmar 
Omega 
Polymax 
Polymax 
Spectrum 
Spectrum 
Suporte 
Sycomig 
Unitron 
Victor do Brasil 
Victor do Brasil 
Microcraft 
Microdigital 
Spectrum 
Gradiente 
Sharp 


MODELO 


Thor 2010 
MC-4000 Exato 
Absolutus 


| Polaris 


DGT-AP 
D-8100 
ENIACII 
Franklin 
Houston AP 
DMII 
MX-2001 
MX-48 

MX.64 
Maxitronic | 
Crafll Plus 
Apple Il Plus 
Apple Master 
Apple Senior 
MC-400 
Maxxi 

Poly Plus 
Microengenho | 
Spectrum ed 
Venus ll 

SICI 

API 

Elppa ll Plus 
Elppa Jr. 
Craftlle 
TK-3000 He 


"Microengenho ll 


Expert GPC4 
Hotbit HB-8000 


FABRICANTE 


Appletronica 
Apply 

CCE 

CPA 

CPA 
Codimex 
Digitus 
Digitus 
Digitus 
Dismac 
Dismac 
Dismac 
Dynacom 
ENIAC 
Engebras 
Filcres 
Pranklin 
Gradiente 
Houston 
Kemitron 
LNW 

LZ 

Magnex 
Maxitronica 
Maxitronica 
Maxitronica 
Maxitronica 
Microcraft 
Microcraft 
Microdigital 
Microdigital 
Microdigital 
Microdigital 


MODELO 


Thor 2010 
Apply 300 
MC-4000 Exato 
Absolutus 
Polaris 
CS-6508 
DGT100 
DGT-1000 
DGT.AP 
D-8000 
D-8001/2 
D-8100 
MX1600 
ENIACII 
AS+1000 
NEZ-8000 
Franklin 
Expert GPC1 
Houston AP 
Naja 800 
LNW-B0 
Color 64 
DMI 
MX-2001 
MX-48 
MX-64 
Maxitronic | 
Craft ll Plus 
Caftlle 
TK-3000 Ile 
TK-82C 
TK-83 

TK-85 


PAÍS 


Brasil 
Brasil 
Brasil 
Brasil 
Brasil 
Brasil 
Brasil 
Brasil 
Brasil 
Brasil 


Brasil 


Brasil 
Brasil 
Brasil 
Brasil 
Brasil 
USA 

Brasil 
Brasil 
Brasil 
USA 

Brasil 
Brasil 
Brasil 
Brasil 
Brasil 
Brasil 
Brasil 
Brasil 
Brasil 
Brasil 
Brasil 
Brasil 


LINHA 


Apple Il+ 


o Sinclair ZX-81 


Apple ll + 
Apple + 
Apple ll + 
TRS-Color 
TRS-80 Mod.lil 
TRS-80 Mod.lil 
Apple + 
TRS-80 Mod.! 
TRS-80 Mod. 
Apple ll + 
TRS-Color 
Apple ll+ 
Sinclair ZX-81 
Sinclair ZX-81 
Apple ll + 
MSX 

Apple Il + 
TRS-80 Mod.lll 
TRS-80 Mod. | 
TRS-Color 
Apple Il + 
Apple Il + 
Apple Il + 
Apple Il + 
Apple ll + 
Apple Il+ 
Apple lle 
Apple lle 
Sinclair ZX-81 
Sinclair ZX-81 
Sinclair ZX-81 


Microdigital TK-90X Brasil 
Microdigital TKS-800 Brasil 


Sinclair Spectrum | 
TRS-Color 


Microdigital TK-90X 
Timex 2000 


Sinclair Spectrum 
Sinclair Spectrum Timex 


Sinclair ZX-81 
Sinclair ZX-81 
Sinclair ZX-81 
Sinclair ZX-81 
Sinclair ZX-81 
Sinclair ZX-81 
Sinclair ZX-81 
Sinclair ZX-81 
Sinclair ZX-81 
Sinclair ZX-81 
TRS-80 Mod.I 
TRS-80 Mod. | 
TRS-80 Mod. | 
TRS-80 Mod.! 
TRS-80 Mod.lil 
TRS-BO Mod.il 
TRS-80 Mod.lil 
TRS-80 Mod.lil 
TRS-80 Mod.lil 
TRS-80 Mod.lll 
TRS-80 Mod.lil 
TRS-80 Mod.IV 
TRS-80 Mod.IlV 
TRS-Color 
TRS-Color 
TRS-Color 
TRS-Color 
TRS-Color 


Apply 
Engebras 


“Filcres 


Microdigital 
Microdigital 
Microdigital 
Prologica 
Ritas 

Timex 
Timex 
Dismac 
Dismac 
LNW 

Video Genie 
Digitus 
Digitus 
Kemitron 
Prologica 
Prologica 
Sysdata 
Sysdata 
Multix 
Sysdata 
Codimex 
Dynacom 
LZ 
Microdigital 
Prologica 


Apply 300 
AS41000 
NEZ-8000 
TK-B2C 
TK-83. 
TK-85 
CP-200 
Ringo R-470 
Timex 1000 
Timex 1500 
D-8000 | 
D-8001/2 
LNW-80 
Video Genie | 
DGT100 
DGT-1000 
Naja 800 
CP-300 
CP-500 
Sysdata Ill 
Sysdata Jr. 
MX-Compacto 
Sysdata IV 
CS-6508 
MX-1600 
Color 64 
TKS-800 
CP-400 


Milmar 
Milmar 


“Milmar 


Multix 
Omega 
Polymax 
Polymax 
Prologica 
Prologica 
Prologica 
Prologica 
Ritas 

Sharp 
Spectrum 
Spectrum 
Spectrum 
Suporte 
Sycomig 
Sysdata 
Sysdata 
Sysdata 
Timex 
Timex 
Timex 
Unitron 
Victor do Brasil 
Victor do Brasil 
Video Genie 


Apple ll Plus 
Apple Master 
Apple Senior 
MX-Compacto 
MC-400 

Maxxi 

Poly Plus 
CP-200 

CP-300 

CP-400 

CP.500 

Ringo R-470 
Hotbit HB-8000 
Microengenho | 
Microengenho Il 
Spectrum ed 
Venus ll 

SICI 

Sysdata Ill 
Sysdata IV 
Sysdata Jr. 
Timex 1000 
Timex 1500 
Timex 2000 
APII 

Elppa Il Plus 
Elppa Jr. 

Video Genie! 


Brasil 
Brasil 
Brasil 
Brasil 
Brasil 
Brasil 
Brasil 
Brasil 
Brasil 
Brasil 
Brasil 
Brasil 
Brasil 
Brasil 
Brasil 
Brasil 
Brasil 
Brasil 
Brasil 
Brasil 
Brasil 
USA 

USA 

USA 

Brasil 
Brasil 
Brasil 
USA 


Apple Il + 
Apple ll + 
Apple Il + 
TRS-B0 Mod.IV 
Apple Il + 
Apple Il + 
Apple l+ 
Sinclair ZX-81 
TRS-80 Mod.Ill 
TRS-Color 
TRS-80 Mod.lll 
Sinclair ZX-81 
MSX 

Apple Il+ 
Apple Ile 
Apple l+ 
Apple Il + 
Apple Il + 
TRS-80 Mod. 
TRS-80 Mod.IV 
TRS-80 Mod.lll 
Sinclair ZX-81 
Sinclair ZX-81 


Sinclair Spectrum 


Apple Il + 
Apple Il + 
Apple ll+ 
TRS-80 Mod. | 


EN an Tin macas 
é, UM LOGUTIE SI qi , 


“INPUT foi especialmente projetado para PN 
“ microcomputadores compatíveis com as sete principais Sinclair E TARSO TK:2000 MSX 


“linhas existentes no mercado. 
ba, então emp a sr 


“Os blocos de textos e listagens de programas po pera 
Eme renes E UM [06] Semp e 
; identificados por meio dos seguintes sim os: TRS-Color pune dl 


ED Dm me eee mo meto mm o em em 





NHRnEnNO PRÓXIMO NÚMERONHEHEEII 


PROGRAMAÇÃO BASIC 


Simulação do balanço de uma empresa. Variáveis 
de distribuição normal. Mercado. Controle de caixa. 


PERIFÉRICOS 


Funcionamento do videotexto. Informações e serviços. 
Conexão ao microcomputador. Compatibilidade. 


PROGRAMAÇÃO BASIC 
Desenhos em três dimensões. Construção do perfil. 
Como alterar o ângulo de visão. Rotação da figura. 


CÓDIGO DE MÁQUINA 


Avalanche: acrescente ao programa a rotina 
que faz Willie saltar à frente. 





